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We propose a nonlinear Schrodinger equation in a Hilbert space enlarged with an ancilla such that the partial 
trace of its solution obeys to the Lindblad equation of an open quantum system. The dynamics involved by this 
nonlinear Schrodinger equation constitutes then a purification of the Lindbladian dynamics. This nonlinear 
equation is compared with other Schrodinger like equations appearing in the theory of open systems. We 
study the (non adiabatic) geometric phases involved by this purification and show that our theory unifies 
several definitions of geometric phases for open systems which have been previously proposed. We study the 
geometry involved by this purification and show that it is a complicated geometric structure related to an 
higher gauge theory, i.e. a categorical bibundle with a connective structure. 
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I. INTRODUCTION 


Quantum information^ and open quantum systems^ are subjects of particular interest in the modern physics, 
dealing with decoherence processes, quantum computation and communication, entanglement processes, distillation 
protocols, Schmidt decomposition, Markovian and non-Markovian effects, etc. A particular interesting subject in 
quantum information is the process of purification^ which consists for a mixed state p of the Hilbert space TLs to 
find a pure state U/ € Hs <8> Ha in an enlarged Hilbert space such that p = tr-^^ | ))«|. The auxiliary Hilbert 
space TLa can be viewed as describing an effective environment. In this paper, we want to study the purification 
process with respect to the dynamics of mixed states. When the dynamics is conservative, i.e. when it is described 
by a Liouville-von Neumann equation ihp = [Hg,p\ (no relaxation effect occurs), the dynamics of the purification 
and the related mathematical structures have been extensively studied, see for example ref>^£. We study in this 
paper the case where the environment of the quantum system induces relaxation effects, with a the dynamics obeying 
to a Lindblad equation^. We will show that the purified state obeys in the enlarged Hilbert space to a nonlinear 
Schrodinger equation. The emergence of a nonlinearity is not a new phenomenon in the relation between dynamics of 
mixed and pure states. In ref£ it is shown that the pure dynamics closest to the Lindblad dynamics (in the sense that 
this pure dynamics is viewed as the dynamics of some tangent vectors on the density matrix manifold) is nonlinear; 
and in ref<£, a purification protocol needed a nonlinear operation has been proposed (a purification protocol is a set 
of operations and measurements transforming a mixed state p to a pure state ip of the same Hilbert space without 
the trace operation in fact the state ip depends only on the protocol and is independent from the initial mixed state 
p -, it is a question different from the purification process discussed in the present paper). Moreover the Liouville 
equation for a piecewise deterministic process is associated for its deterministic part with a nonlinear Schrodinger 
equation but which does not take into account the jump part (see ref^ chapter 6.1) in contrast with our equation for 
the purified dynamics. 

Geometrization of physical theories is a great active area in theoretical physics. In nonrelativistic quantum dynam¬ 
ics, it is in particular related to the theory of geometric phases (so-called Berry phases)^. As shown by Simon^, 
the dynamics of pure states and the Berry phases take place in a geometric structure which is a principal fibre 
bundle endowed with a connection. Some generalisations of geometric phases have been proposed for mixed states: 
UhlmanniSr— has proposed a concept of geometric phases based on the theory of transition propabilities for pair of 
mixed states, the involved geometric structures have been analysed in ref— ; Sjoqvist eta hf&^ have proposed a 
concept of geometric phase based on an interferometric theory, the involved geometric structures have been analysed 
in reGid£; and we have proposed a concept of geometric phase in the adiabatic limit based on the generalisation 
of the geometric structure studied by Simon from vector bundles to C^-module s ^ ' —. By using the equation of the 
purified dynamics, we will built a general theory of geometric phases for open quantum systems which unifies these 
previous approaches. The dynamics of the density matrices present different geometric phases appearing at different 
levels. This is due to a more complicated gauge structure called higher gauge theory in the litterature2I~—. We will 
show that this struture is a generalization in the category theory of the principal bundle structure, complicated by 
the stratified structure of the density matrix manifold^. 

This paper is organized as follows. Section II is devoted to the geometry associated with the purification process 
by recalling the stratified structure of the density matrix manifold and by introducing some representations of the 
purified states with the associated inner products. This section introduces some mathematical tools needed to the 
understanding of the following. Section III shows that the purified state satisfies a nonlinear Schrodinger equation 
if the associated mixed state satisfies a Lindblad equation. A general theory of geometric phases for open quantum 
systems is discussed, and there relations with the different propositions of geometric phases are analysed. An example 
consisting to a two level system driven by an external field and submitted to dephasing and spontaneous emission 
effects is discussed with a numerical integration of the nonlinear Schrodinger equation and numerical calculations 
of geometric phases. Section IV studies the geometric structure involved by the general geometric phase theory 
and the purification process, firstly from the point of view of the ordinary differential geometry, secondly from the 
point of view of the category theory. It concludes by the introduction of the connective structure and the physical 
interpretations of the different fields involved by the connection. 

A note about the notations used here: 

We adopt the Einstein’s convention: a bottom-top repetition of an index induces a summation. 

B{TL) denotes the set of the bounded linear operators of the Hilbert space H. For an operator A £ B(TL), RanA, ker A 
and Sp(A) denotes its range, kernel and spectrum. VA, B E B[fH), we denote the commutator and the anticommutator 
by [A,R] = AB — BA and {A,B} = AB + BA. Autg, with g a vector space or an algebra, denotes the set of the 
automorphisms of g. 

The symbol between two spaces (resp. manifolds) denotes that they are isomorphic (resp. homeomorphic). The 
symbol “ss ” between two manifolds denotes that they are locally homeomorphic. The symbol ‘A B ” denotes an 
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inclusion of A into B. G xi H denotes a semi-direct product between two groups G and H; g 5 f) denotes a semi-direct 
sum between two algebras g and t). 

Pr^ : V\ x V 2 x ... x V n —> V., with Vj some sets, denotes the canonical projection Pr^i, X 2 , ■■■,x n ) = Xi. Let M be 
a manifold, T X M denotes its tangent space at x (TM denotes its tangent bundle) and Ll n (M,g) denotes its space of 
g-valued differential n-forms. Let f : M —> N be a dijfeomorphism between two manifolds, /* : TM TN denotes 


its tangent map (its push-forward) and f* : Ll*N —> LI*M denotes its cotangent map (its pull-back). Let E M be 
a fibre bundle (E and M are manifolds and tv is a surjective map); T(M,E) denotes the set of its local sections. 

For a category ^, Obj'tf denotes its collection of objects and Morph‘d denotes its collection of arrows (morphisms). 
Vo £ Obj^, id 0 denotes the trivial arrow from and to o (the identity map of o). Va £ Morph^, s(a) denotes the 
source of a, andt(a) denotes the target of a. Funct (‘if,'^') denotes the set of functors from toff' fEndFunct(‘^ ) ) = 
Funct(^, *£)). 


\/A(t) £ B(H), Te~ f ^ A{t ' )dr 


Ua( t, to) denotes the time-ordered exponential (the Dyson series) i.e. the solution 


of the equation: dUA ^’ to ) = —A(t)UA(t,to) (with UA{to,to) = id n)- Te ^0 A( ' t ' ,dt = denotes the time-anti- 

ordered exponential, i.e. the solution of the equation: dVA< d t °’ t ' > = —VA{to,t)A(t) (with VA(fo,to) = vA-h)- 


II. PURIFICATION PROCESS 

A. The purification bundle and the stratification 

Let Hs be the Hilbert space of the studied system S (for the sake of simplicity we suppose that Hs is finite 
dimensional, Hs — C"). Mixed states of S are represented by density operators; the space of density operators being 

L> 0 (H s ) = {p £ = P, P > 0, trp = 1} (1) 

where B(Hs) — 9Jt„xn(C) denotes the space of linear operators of Hs- The mixed state oj p : B(Hs) —> C associated 
with p £ T>o(Hs) is defined by 

VA € B(H S ), oj p (A) = tr(pA) (2) 

B(H S ) being a C*-algebra and w p being a state of this C*-algebra (see refi^). In some cases, we can also consider 
non-normalized density operators: 


T>(H s )={p£B{H s ),p'=p,p> 0} (3) 

VA £ B{H S ), w p (A) = ^M (4) 

trp 

A purification of a density operator p £ T>(Hs) is a state 'Lp £ Hs <8> Ha such that 

p = tr w J'F p »«'I , p | (5) 

where Ha is an arbitrary auxiliary Hilbert space called the ancilla and ((.|.)) denotes the inner product of Hs < 8 > Ha 
induced by the tensor product (tr-^ denotes the partial trace over Ha)- Because of the Schmidt theorem^ one needs 
dim'H.A > dim'Hs. In order to avoid any unnecessary complication, we choose AvovHa = dim'Hs throughout this 
paper. 

An interesting choice of ancilla consists to the algebraic dual of Hs, Ha = Hs* = € Hs} (the space of 

continuous linear functionals of Hs). In that case Hs < 8 > Ha = Hs ® Hs* = B(Hs)- This choice is called standard 
purification. B(Hs) is endowed with the Hilbert-Schmidt inner product: 

VW,Z £ B(H S ), (Z\W) HS = tx Hs {Z'W) = Z ta W ia ( 6 ) 

where W = W ,0 j£.j)(C a |, (Ci)i=i,...,n being an orthonormal basis of Hs- We have then 


tvhs ■ B(H s ) -f D(H s ) 

W 1 — ww f 


( 7 ) 
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FIG. 1. Simplex E(n) for n = 2 (left), n = 3 (middle) and n = 4 (right). The strata are the open sub-simplices of each simplex 
(including E°(ra) the interior of E(n) - the stratum of no degenerated eigenvalues -). For each stratum E 1 the model manifold 
of 7r)) 1 (a') is indicated. Singular strata are written in grey whereas regular strata are written in black (• denotes the stratum 
of the microcanonical state, and o denotes the stratum of the pure states). The description of the simplexes can be found 
table [H The details about the stratification can found in ref.' 3 where E(n) is called Weyl chamber by analogy with the group 
representation theory. 


The restriction of 7 ths on ShsB{Hs) = {W £ B(Hs), ||W||^ S = tv-H s {W^W) = 1} has its values in T> 0 (Hs). n hs is 
a surjective map. 

MW €S hs B(Hs),VU €U(U s ), W = WU=>tths(W) = WUU*W* =k hs {W) (8) 

U(Hs) — U(n ) denotes the group of unitary operators of Hs- The action of U(Hs) on ShsB(HLs ) is not transitive 
except for the restriction on the faithfull operators (det W ^ 0). 

U[HLs) ac ts also on the left of Shs^CB-s)- This action induces the adjoint action of the group on V 0 (T-Ls ): 

MW £ S hs B(H s ),VU £ U(H S ), W =UW ^ir HS {W) = Uir H s(W)U~ 1 (9) 

Vp £ 'DoCHs), the orbit of the density operator U('Hs)p = {UpU -1 , U £ (7(Ps)} is constituted by all operators which 
are isospectral to p. T,(Hs) = 'D 0 ('Hs)/U('Hs ) — £(n) where 

n 

T,(n) ={{pi,...,p n ) £ [0,1 ]",Pi <Pi+i,^Pi = 1} (10) 

i=i 

is the (n — l)-simplex of the possible eigenvalues (dim® £(n) = n — 1). Y>(Jis) is the set of diagonal density operators 
with sorted eigenvalues. Let 7 td ■ Do^H-s) ^{HLs) be the canonical projection associated with the quotient space 
V 0 (H s )/U(Hs). 


Va e E(H 5 ), t to 1 (a) = {UaU~\U £ U{U S )} ~ U{ns)/U{n s ) a (11) 

where U(Hs)a = {U £ U(Hs), UaU ~ 1 = a} is the isotropy subgroup (the stabilizer) of a. The bundle V 0 (Hs) 
£(7£s) is not locally trivial, it is a stratified space^^i. A stratum E ! (Ks) of £( H-s ) is characterized by the degeneracy 
profil of the eigenvalues. If a £ £ z ( Hs) has k\ eigenvalues with degeneracy equal to qi, /c 2 eigenvalues with degeneracy 
equal to q 2 , etc, then its stabilizer is U('Hs)a — U{q\) kl x ... x U(qi) kl (the products are the group direct products) 
and then vr” 1 ^) ~ U(n) / (U{qi) kl x ... x U(qi) kl ) Mi = n). The stratum T,°(Hs ) with no eigenvalue degeneracy 
is associated with the normalizer U(Hs)(j — T n (the n-torus) and Vcr S E°(Ks), 7rjj 1 (a) ~ U{n)/T n = Fl(n,C) (the 
flag manifold); in the other cases, the fibers of the strata in Do(Hs) are homeomorphic to partial flag manifolds as 
Grassmanian manifolds Gr p ( C") = U{n)/(U{n — p) x U{p )), projective manifolds CP" -1 = U(n)/(U(n— 1) x (7(1)), 
etc. We distinguish between the regular strata such that det p ^ 0 (no zero eigenvalue) and the singular strata such 
that det p = 0. We remark the presence of two particular strata associated with vertices of £(n), the singular stratum 
of the pure states £ °{'Hs) (pi = ■■■ = p n -1 = 0 and p n = 1, n^^a) ~ U(n)/(U(n — 1) x (7(1)) = CP" -1 ), and 
the regular stratum of the microcanonical state £°°(7is) (pi = ... = Pn = T, 7 rjj^cr) — U(n)/U(n) = {1})- Fig. [Q 
illustrates the stratification for the smallest dimensionnal cases. 
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TABLE I. Descriptions of the stratification of the simplex E (n). 


n 

strata E 1 (n) 

sub-simplices 

fiber (a) 

2 

0 < pi < p 2 

edge 

Fl( 2,C) 


CN 

i—l 

II 

<N 

II 

s 

vertex 

{1} 


Pi = 0 < P2 = 1 

vertex 

CP 1 

3 

0 < Pl < p 2 < P3 

facet 

PZ(3, C) 


0 < pi = p 2 < P3 

0 < pi < p 2 = P3 

edges 

CP 2 


Pl = 0 < p 2 < P3 

edge 

Fl{ 3,C) 


Pl = 0 < p 2 = P3 = 1/2 

Pl = P2 = 0 < p 3 = 1 

vertices 

CP 2 


Pl = P2 = P3 = 1/3 

vertex 

{1} 

4 

0 < pi < P2 < P3 < P4 

cell 

PZ(4, C) 


0 < pi = p 2 < P3 < P4 

0 < pi < p 2 = P3 < P4 

0 < pi < p 2 < P3 = P4 

facets 

Z7(4)/(Z7(2) x T 2 ) 


0 = pi < p 2 < P3 < P4 

facet 

PZ(4,C) 


0 < pi = p 2 < P3 = P4 

edge 

Gr 2 ( C 4 ) 


0 < pi = p 2 = P3 < P4 

0 < pi < p 2 = p 3 = P4 

edges 

CP 3 


Pl = 0 < p 2 < P3 = P4 

Pl = 0 < p 2 = P3 < P4 

Pl = P2 = 0 < p 3 < P4 

edges 

Z7(4)/(t/(2) x T 2 ) 


Pl = P2 = P3 = 0 < p 4 = 1 
Pl = 0 < p 2 = P3 = P4 = 1/3 

vertices 

CP 3 


Pl = P2 = 0 < p 3 = P4 = 1/2 

vertex 

Gr 2 (C 4 ) 


Pl = P2 = P3 = P4 = 1/4 

vertex 

{1} 


Let Ha be an arbitrary ancilla, and let (C*)i=i,...,n and (€a)a=i,...,n be orthonormal basis of Hs and Ha- Let 
: B(LLs) — > Lis <8> Li a be the (non canonical) isomorphism defined by 

VIE G B(U S ), *{W) = W ia Ci ® (12) 

with W = lE* a |Ci)(Ca|- The restriction of \l/ on ShsB{Hs) has its values in S(Lis ® Ha) = G Lis ® 
Ha > = 1} ~ 5 ,2n2 - 1 (the (2n 2 — l)-dimensional sphere). Note that ((’® , (Z)|S1/(1E))) = (. Z\W)hs■ 

The following commutative diagramm summarizes the bundle of purification: 

ShsB(H s ) —J-> S(H s ®Ha) — S 2n2 ~ l 

™ s j. ”*{ 1 

ZWs) - 2>o(fts) —2? 0 (C") - Ui®&(C") U(n)/U(n) a 




-i- Np 



X(Hs) - S(Ws) — t E(n) - UiS^n) f — M 

where denotes inclusion maps and ~ denotes homeomorphisms, the vertical arrows are projections. 

The ancilla plays the role of an effective environment with which the system is entangled. The right action of U(Hs) 
on ShsB(Hs) is associated with unitary transformations of Ha- 

MU G U(H S ),\/W G S HS B(Hs), 'SHWU) = l Hs <g> U T 'S>(W) => ^{^{WU)) = p w (13) 

where U T = (C j g|Z7C a }|£ a )(£ /3 | € U(Ha) ( t denotes the transposition) and pw = 7r^(^(lE)) = khs(W) = WW^ . 
Such a transformation has no consequence for the observator which has only information concerning the system 
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(information concerning the environment is lost by the partial trace tr%^). It is then associated with an inobservable 
reconfiguration of the ancilla (the effective environment). In contrast, the left action of U(Hs) on ShsB{Hs) is 
associated with unitary transformations of Hs modifying the mixed state: 

VU€U{Hs),VWeS H sB{H s ), ^(UW) = U®\ UA ^{W)=^^(<b(UW)) = Up w U- 1 (14) 

In quantum information, it can be interesting to consider also SLOCC transformations (Stochastic Local Operations 
and Classical Communication). By virtue of the principle of the quantum open systems, information concerning the 
ancilla (the environment) is lost, the physicist cannot performs SLOCC transformations on Ha- The group of SLOCC 
transformations on Hs is SL(Hs) — SL(n, C) (the group of invertible operators of Hs with determinant equal to 1). 
Since U(Hs) 3 SU(Hs) C SL{Hs) ( SU(Hs ) is the group of unitary operators of Hs with determinant equal to 1 ), 
the group of unitary and SLOCC transformations are SL(Hs) x 17(1) = GL(Hs) (17(1) is the group of phase changes 
and GL{Hs) is the group of invertible operators of Hs)- 

VU e GL{H S ),VW e B(H S ), ^{UW) = U ® 1 Ha V{W) => n*{'S'{UW)) = Up w U ] (15) 

by noting that ShsB(Hs) and V 0 (Hs) are not stable by the right action of SL{Hs) (we must then consider B{Hs) 
and T>(Hs))- We extend ttd on the whole of D(Hs) by V17 G SL(Hs ), Vp € 2?o(Hs), no(UpW) = i td(p)- 


B. The C’-module structure 

The purification space Hs ® Ha can also be viewed as a left Hilbert S('H l s)-C'*-modul82&, the left action of the 


C*-algebra B(Hs) being defined by 

VAeB(Hs),W eHs®n A , = A®1 Ha ^ = ¥ a {ACi)®^ a (16) 

The inner product of the C*-module is defined by 

V*,Q€Hs®H a , =tr*Jtf»«$| GB(H S ) (17) 

It satisfies the usual properties: 

+ (18) 

V'f,$ G H s ® H 4 , <Wt = <*|$>* (19) 

Vtf G Us <8 Ha, <^1^)* > 0 and (T|T)* = 0 ■<=> T = 0 (20) 

The density operator associated with a state of Hs ® Ha is its square norm with respect to the C*-module structure: 

VTgHs®Ha, 11^112 = (^),=p^GV(H s ) (21) 

B{Hs) as the standard purification space is related to the C*-module structure by 

VHg B(H s ),VW G B{H S ), *(AW) = A^{W) (22) 

VW, Z G B{H S ), <®(W)|^(Z)), = ZW f (23) 

The right action of B{Hs) on itself induces a particular class of operators of the C*-module acting only on Ha- 

VH G B{H s )yW G B{H S ), *(WA) = l Hs ® A T *(W) (24) 


C. C*-adjointness 

In this part we want to define the adjoint, with respect to the inner product of the C*-module, of the right action 

of B{H S ) on H s ® Ha- 
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Definition 1 Let T G B(Hs)- We define Tt : Lis 0 Ha — t B(Hs) by 

v* e h s ®h a , ri(t) = ( 25 ) 

= mrw*) 0 a \t a ){^\ (26) 

= mrw*r a ^ a )(^\ ( 27 ) 

where Wy G B(Hs) is such that ^(W 7 ^) = 5', T l j = (CITCj), \IP a = W^ a = (C^W^C 12 ) = ((0 7 ® C Q |^))> and where 
W G B(Hs) is the pseudo-inverse ofW^: 

IFJ IFj/ = 1 — PkerWy and W$WJ = PRanW* (28) 

where Pk er w* and PRanW* are respectively the projections on the kernel and on the range of Wy in Lis ■ 

Property 1 Let rt(tf) = r + (\I/)^ (the adjoint f is in the sense of the inner product of Li a)- We have 

v$G« s ®^, <$|i„ s ®n(’W* = (i« s ®rt(^)$|'f)* (29) 

«=► (<f>\PTUnw*r®iH^)* = (i'Hs®r t W'S>\*)* ( 30 ) 

For this reason, we call r* the C*-adjoint of T. 


Corollary 1 

Proof: 


V'F G Lis 0 LLa, 1 Hs ® F+ (tp)Vp = PRanW^r ® 1 u A ^ 
1 Hs ® rl(^)^ = T^) a p ^Ci ® 

= r fc J -$ J ' a (wj).pfc^Ci®ea 

= r fc i ^“(W’*w’|) < fc Ci®&« 

= (PRan^) i fe r fc i ^ a C i ®^ 

= Pk&iiW* r 0 l'H s tP 


(31) 

(32) 

(33) 

(34) 

(35) 

(36) 


($|l„ s ® n(tf)tf), = tr„Jl„ s ® H(tt)*»«S| (37) 

= n(^r /3 ^l-|c i )(C j | (38) 

(l ns Ti(*m*)* = tr*JF))«l Ws ®r*($)$| (39) 

= ^ffw 7 i-|Ci)(c j 'i ( 40 ) 

= ^n ( vt)^$-|c i )<c j | (41) 

□ 


It can be interesting to relate RanlF^ (which appears in the definition of the C*-adjoint) with the density operator. 

Property 2 Let ’F G H S ®H A , W* G P(P 5 ) 6e skc/i ¥(W®) = \E>, crndp* = \\^\\l = w A (V) = W^W\ G V(H S ). 
We have RanlhV = Ranp^ and coker IF^ = ker p^ . 

Proof: coker Wp = {</> G Hs, = 0}. <fi G cokerfF^ => 0,'I' 4 “'3/j a = 0 => fifip^) 1 ^ = 0 since py = \IT a \E'j a |Ci)(C' 7 1- 
We have then {<j>\p^ = 0 => py\<f>) = 0 => (f) G ker p#. We have then coker IF^, C kerp^. 

(p*)F|C») = ^“^jalCi) => V0 G Hs, {p^yja^^Ci)- We have then V0 G Hs,p^f> = Wy(Wl<j)) => 
pyj) G RanTF^. We have then Ranp^ C RanlFj,. 

Suppose that ((\)i the basis of Hs is the eigenbasis of p^. (p^)* ■ = = p l 5j. We can write (’Fjl'F 4 )^ = p l 5 *• 

with F' = \IT“£ a g Pkt- We have then p l = It follows that p 1 = 0 => F = 0 => F“ = 0,Va. RanW 7 ^ = 

{F Q 0 Q |£i}, (f> a G C} and then RanlFj C coRanp^ (coRanp^ = (kerp^)- 1 ). But coRanp^, = Ranpjj, = Ranp^ 
(pl = py). We have then RanlFj C Ranp^. 

The two last paragraphs show that RanlFj = Ranp^,. The first one shows that cokerC ker p^, but dim cokerfF^ = 

_L 

dim Hs — dim RanW 7 ^ (by definition Hs = RanW® 0 coker). We have then dim coker= dim Hs — dim Ranp^ = 
dim ker p^ (by virtue of the rank-nullity theorem). This induces that coker Wy = ker p^. □ 



III. PURIFICATION OF LINDBLAD DYNAMICS 


A. The nonlinear Schrodinger equation of the purified dynamics 
1. From the Lindblad equation to the nonlinear Schrodinger equation 

Theorem 1 Let 1 i—X p(t) G T>o(Hs) be a solution of the Lindblad equation 

Ihp = [H s , p] - ]p k {T[r k ,p} + i^ k T k pT\ (42) 

with Hs G B(Hs) the system Hamiltonian, T k G B(TLs) the quantum jump operators and q fc G [0,1] the relaxation 
rates. Let 1 i—X 4/ p (f) G 'Hs'Si'Ha be a purification state of p(t). ^f p is solution of the following projected non-hermitian 


nonlinear Schrodinger equation: 

(1 Hs ® -PRanp)*ft’i'p = (1 Hs ® -pRanp) {lis 0 1 HjC^ P - 0 1 U A ^ p + k ® 0^p)^b>) (43) 

Proof: Let W p (t) G BifHs) be a standard purification of p, p = W p Wj ,, such that SI/(W P ) = Sl/ p . 

ihp = C(p) (44) 

=► tnw p wl + ihw p wl = [ h SiP ]-!7 fc {r!.r fc , P } + 2x | 7 fc r fe prt (45) 

=► thw p wl - (z^w p w'+) t = - ^7 fe rir fcP + l -j k r k pr{ - (h sP - l 7 *rtr fcP + l -i k T kP T \) 1 (46) 

=► *^p^ = \p + |7 fc r fc prt + K (47) 

where K = K' is an arbitrary self-adjoint operator. We can set I\ = 0 without loss of generality. 

ihw P wl = H s w p wl - ^rir.w.w} + ^ k nw p w}r{ (48) 

=► ihw p p coRanWp = H s w p p coRanWp - ^ k rlT k w p p coRanWp + l -j k T k w p wlv\{wlY ( 49 ) 

By application of SI/ on this last equation we find 

1 Hs ® PlUnplMp = H S 0 XRanp^p - ^T^k 0 P R anp’J'p + ^k 0 (^^ (Wj )*) T X p (50) 

But (wt r t ( wt)*) T = (w;iw p )t T = r|(x p ). □ 

On the regular strata, no projection occurs and the purification state is solution of the nonlinear Schrodinger 
equation: 

im p = (h s - \l k T\T k ) 0 1 u^p + \l k T k 0 r|(X p )X p (51) 


2. From the nonlinear Schrodinger equation to the Lindblad equation 

Theorem 2 Let SH G TLs 0 BLa be a solution of the nonlinear Schrodinger equation: 

iMt = (h s - k) 0 1 + ^7 fc r fc 0 r|(tf)tt (52) 

Then p ^ = 7T^i(X) = ||X||^ G V(TLs) Is solution of the following master equation: 

ihpq, = [H s ,p*] - ^7 fe {r^r fc ,p^} + ^7 fe {r fc p^r[,p Ranp ^} 


( 53 ) 



9 


Proof: Let H u {.) = (H s - |7 fe rl r fc) 0 1 u A + f 70 r|(*) be the nonlinear operator of the purified dynamics and 

/V- 'I'M*. 

ihPn, = | (54) 

= \H u (*)n({*\ - (55) 

= [Hs ® 1 h a ,p*] - ^ 7 fe {r[.r fe ® 1 h a ,p*} + ^ 7 fc r fe ® r|(*)p* + | 7 fc ^r fc 0 r* (*) (56) 

By applying the partial trace tr-^ on this last equation, we find 

ihp* = [Hs,pv\ - \i k {T\Tk,p*} + ^ 7 fc (^|r fe 0 r* + ^ k (r k 0 r|(*)*|*>. ( 57 ) 

Since T k 0r|(\H) = (1 n s 0 r|($))(r fc 0 1 n A ) we have 

*N>* = lH Sl p *] - \i k {r\T k ,P^} + \i k {iu s ® r£(*)*|r fe 0 i n ^)* + ^7 fc (r fc 0 i Ha ^\Ih s ® r* (*)*>. ( 58 ) 

= [Hs,Pv] - t^IT^T*,,^} + -j k (PRanp^k 8 8 1 Ua^)* 

+ 2^ k ^ k ® Iw^^l-PRanp^rfe 0 lu A ^)* (59) 

= [Hs,p*] - l -i k {T\r k , p*} + iy^r^rtiWp, + ^7 fe ^Ra„p*r fcP>P r[ (60) 

□ 

On the regular strata, i-Ranp* = 1 -h s and the master equation reduces to the Lindblad equation: 

ihpM, = [H s ,p<z] - ^7 fc { r l,r fc ,p^} +i'y k T k p^,v\ (61) 


3. Comparison with other Schrodinger equations associated with the open quantum systems 

We have shown that the purified state obeys on the regular strata to the nonlinear Schrodinger equation: 

ih*p = HHf®l n ^ p + ^ k r k ®Tl(* p )V p * p eHs®H A (62) 

with H e H = Hg — |r),rfe. Reb£ studies a nonlinear Schrodinger equation for the pure evolution closest to the 
Lindblad evolution in the sense that it is viewed has the dynamics of some tangent vectors of T>o(Hs)- This equation 
is the following: 

iM = (H e H - (H e ff)^ + * 7 fc (r[)p(r fc - (r fc )p)^> 4 e u s (63) 

where (0)$ = {ip\0\ip} (VO € B(Hs))- Except the shift of the operators by their average values, the structure of this 
equation is similar to eq. (16211 . In fact suppose that during a short time the dynamics remains on the singular stratum 
of the pure states: p = \ip)(ip\, in that case W p = \ip){ipo\ and W* = \ipo)(ip\ (with ipo = ip(t = 0)). The equation (l43l) 
becomes then 

ihip 0 '0 = (7L e// 0 + ^7 fc ( r fc)^r fc 0) 0 ip (64) 

which is very close to the equation (lixil) (except the shift of the operators and a | factor). 

The deterministic part of the Liouville equation for a piecewise deterministic process is (see rei£ chapter 6.1) ihip = 
H e Hip -f | 7 fc ||r fc 0|| 2 0- Since ||Tfc|| 2 = (r^Tfc)^,, it is very close to our equation for a dynamics remaining on the 
pure state stratum. This is valid only between two quantum jumps; in a same manner equation (1641) could be 
approximatively valid only during a short time until the quantum jumps described by the operator lp'y k (r]f) k tend 
to leave the pure state stratum. To describe quantum jumps we add a stochastic process (see refi2>2£): 

ihdip = H e Hipdt + ^ k (r{Y k )^dt + idN k - lw*J 4 4’£Hs (65) 

or ihdip = H eff ipdt + ^ 7 fe ^(rj. + T k )^T k - ^(rj. + T fe ) 2 ^ ipdt + ly/j^dWf ^T fe - i(T|. + T k )^\ ip (66) 
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TABLE II. Properties of the Schrodinger equations associated with an open quantum system. 


Schrodinger equation 

Generator of the quantum jumps 

Nature 

dim. 

relation with the density matrix 

Hilbert-Schmidt representation 

fc ® rf 

“superlinear”^ 

n 2 

H 

II 

Q. 

Purification 

^7 fe r fc ®rt(^ P ) 

nonlinear 

n 2 

p = trwJ'i'pM'i'pl 

Closest pure evolution 


nonlinear 

n 

tr ( Sp — 5\ip){il>\) 2 is minimum^ 

Piecewise deterministic process 

".‘(lirA, 

nonlinear stochastic 

n 

P = E (WW) 

Quantum state diffusion 

i^dw?{r k -\{r k + rl)^) 

nonlinear stochastic 

n 

P = E (WW) 


a linear in the sense of the superoperator theory 


where dN k is a Poisson process and dW k is a Wiener process. The first equation is the representation as a piecewise 
deterministic process and the second one is called quantum state diffusion process. The density matrix is then 
p = E (|'0)('i/>|) where E is the expectation value with respect to the stochastic process. 

Ref . 30 proposes a Schrodinger equation following the dynamics induced by the Lindblad equation: 


ihT p = H e fi ®lu^ P -lus®{H e ff)^r p +i 1 k T k ®T^r p r p €H s ®H A (67) 


(the operators of Hs being transformed into operators of Ha by the non canonical isomorphism O a b|Ca)(C b | h A 
O a p\£,u){£,^\ ((C a)a and {£ a )a being the two choosen fixed basis of Hs and Ha))- Equation (1671) is linear, but it is 
not an equation for the purification, since in this case T p = p a pCa ® ^ and then tr?^ | Y p ))((Y p | ^ p. In fact, if we 
choose the Ha = Hs* , equation (1671) is precisely the Lindblad equation in the Hilbert-Schmidt representation (in 
the Liouville space), |T p )} = p a b \( a ) ® (C 6 | = P (equation (l67l) is then just a reformulation of the “superoperator” 
formalism). 

The comparisons between the different Schrodinger equations associated with an open quantum system are summarized 
table HIl 


4. Example: control of a two-level system with dephasing and spontaneous emission 

a. The model: We consider the two level system governed by the following Lindblad equation 

ihp= [Hs , p\ — ]p k {T\T k: p}+ i^ k T k pT\ ( 68 ) 

with 

H s (t) = -pS-B(t) (69) 

T *(t) = R(B(t))a z R(B(t))' (70) 

T_(f) = R{B(t))a-R{B{t))^ (71) 


where S = -|(j, <r_ = a x — i<j y = ^ ^ ^ J, {a X: a y , a z } are the Pauli matrices, and R(B) is the matrix of the 

eigenvectors of Hs expressed in the canonical basis of C 2 (the effects of the environment must be described on 
the dressed basis of the system, see for example refi^L). This model corresponds to a spin-1/2 interacting with 
a magnetic field B [p being the coupling constant) with an environment inducing a dephasing modelled by r 2 (t) 
and a spontaneous emission from (°) to (°). 'y z and 7 “ are the square of the probabilities per unit time of these 
decoherence processes. Hs = C 2 and then Hs®Ha = C 4 . We can rewrite the operators by using spherical coordinates 
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FIG. 2. Strenght B and colatitude 8 of the magnetic field controlling the spin, with respect to the time t (atomic units). The 
longitude ip is remained to zero. 


(.B = B (sin 8 cos <pe x + sin 8 sin ipe y + cos 8e z )): 


H S = - 


R = 


r z = 


r_ = 


fihB I cos 0 e lv sin 8 
2 \ e lv> sin 8 — cos 9 


cos | 
e ltp sin 




e 


cos 


e 


2 2 

cos 8 e~ llp sin 8 
e lLp sin 8 — cos 8 


-e-^sin 8 -2e" 2 ^ sin 2 § 
2 cos 2 | 


e l(p sin 8 


(72) 

(73) 

(74) 

(75) 


b. The dynamics: We consider the evolution of the magnetic field controlling the spin, which is represented figure 
[2] We consider the initial condition 


P( 0) = 


6 e 
6 1 — e 


(76) 


with e <C 1 a very small parameter. We choose this initial condition in order to be close to (” J) but without starting 
from the singular stratum of the pure states 2?°(C 2 ) (we avoid then the problems induced by the transition from a 
stratum to another one). The dynamics runs in the regular stratum T>°( C 2 ), but the steady state induced by the 
Lindblad equation is in the singular stratum 2?°(C 2 ): 


lim p(t) 

t->-+oo 



(77) 


because of the spontaneous emission process. The evolution of p(t) solution of the Lindblad equation and the evolution 
of 4+) solution of the nonlinear Schrodinger equation = {{Hs — | 7 ^c ^] c ^fe)®lc 2 + | 7 fe ^fe®^|('^'))\^ , are represented 
figure[3] We can note that p(i) = 7r+\I/(t)) as shown hgure|4j The dynamics shown figure[3]can be compared with 
the other Schrodinger dynamics associated with the system shown figure [5] 


B. Geometric phases of open quantum systems 
1. The different notions of operator valued phases 

For the open quantum systems, the notion of phase cannot be the same that for the closed systems. In accordance 
with the C*-module structure of the purification space TLs <8> Ba, a phase for an open quantum system is certainly 
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FIG. 3. Left: evolution of the populations and of the coherence of the mixed states p(t) (solution of the Lindblad equation). 
Right: evolution of the occupation probabilities of the purified state 4/(t) (solution of the nonlinear Schodinger equation). 
(l*J l }))ib=i ,2 is the canonical basis of Hs ® Ha = C 2 ® C 2 . The parameters are chosen as being (in atomic units) p = 1, 
-y z = 10 — 3 , 7 ~ = 2 x 10~ 3 , and t = 10~ 5 . The Lindblad equation is integrated by using a splitting of the evolution operator 
U(t, 0) = Te - ^ Jo £(t )<u —jj(j- N , t N _i)...U(t 2 ,ti)U(ti, 0) (for a partition 0 < ti < ... < tjv of [0, T]) with the approximation 

U(ti+i,ti) = Te~ th b c W dt — e - thC ( t i') At + 0(At 2 )-, where £ = — 37 *^]^, •} + * 7 *T fc »r] s is considered as a linear 


(non-hermitian) operator of the Hilbert-Schmidt space Hs ®Hs* (as a “superoperator” acting on the Liouville representation 

of p : [ ei 2 j ). The nonlinear Schrodinger equation is integrated by using an intermediate representation 'F(f) = U e ^{t, 0)~ J ® 
V P22 J 

1 C 2 'F(t) with = [Hs—^ k V\Tk)U e ^ integrated by using a splitting on the time partition, and by using a RK4 (Runge- 

Kutta of order 4) algorithm to integrate the nonlinear part Jj'F = ^-•y k U e ^~ 1 YkU e ^ <g> ^(U e ^W^)Wl.(U e ^W ^,)*^j 
with Wy = <I’~ 1 ('I'). 



FIG. 4. || p(t) — 7r^(4/(t))|| 00 in logarithmic scale (with ||M||oo = maxij \Mij\), p(t) being obtained by integration of the 

Lindblad equation by a split operator method and ^(t) being obtained by integration of the nonlinear Schrodinger equation 
by a mixing of a split operator method and a RK4 algorithm (as explained in the caption of figure [3]). Note that the step of 
the time partition is At = 0.225 (atomic unit). 


an operator in the C*-algebra B(Hs) (we recall that a C*-module mimics the structure of vector space by replacing 
the ring C by a C*-algebra). 

Definition 2 We define two different notions of phase in the C*-module for a state U/ £ Hs <§) Ha' 

• We call phase by invariance, an operator g £ GL(Hs) leaving invariant the C*-norm: 

115*112 = 11*11* (78) 

• We call phase by (unitary) equivariance, an operator g £ GL(Hs) leaving equivariant the C*-norm: 

bn^gimlg- 1 


(79) 
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FIG. 5. Up: Solutions of the nonlinear Schrodinger equations of the closest pure dynamics with lhtp pu r = (H e ^ — 
(H eff )^ pur )ip P ur + I 7 fc (rj,)y, pur .(r k - (T k )^ pur )ip pur (left) and ihtp pur = H eff ip pv , r + §7 fc (Tl)p pur Tkippv.r (right). Up to 
a renormalisation factor, these evolutions are very close to the evolution of the isolated system ( 7 * = 0). Down: So¬ 
lutions of the stochastic Schrodinger equations for the piecewise deterministic process with ihdippDP = H e ^ipPDpdt + 
h 1 k {^k^k)Tijp D pipPDpdt + dNt ( ||r fc i^p DP || — ipPDP (left) and for the quantum state diffusion with ihdipQSD = 

H*”j, QgD dt+Z 7 fc (<rt + r^ QSD r fc - £<r£ + r fe )^ QSD ) ^Qsndt + uffidW? (r fc - §<rj + r fc )^ QSD ) <Pqsd (right), E de¬ 
notes the expectation value with respect to the stochastic process. The stochastic dynamics are computed with a Monte-Carlo 
method with 30 trajectories. 


For an abelian C*-algebra (as C) the two notions are equivalent. We can note that Vg £ GL(Jis ), V’F £ Hs we 

have 


110 * 112 = 511 * 11 ^ ( 80 ) 

and then each (SLOCC or isospectral) tranformation of S appears as a phase by non-unitary equivariance (a concept 
similar with the case of the non-hermitian quantum systems: due to the non-conservation of the norm during the 
dynamics, we consider “non-unitary phases” g £ C* corresping to a change of norm:||cp/>|| 2 = |</| 2 ||'0|| 2 with if) £ Tis)- 
Let p £ V{7is)- The group of phases by invariance of 7r) 4 1 (p) is the stabilizer (the isotropy subgroup) of p for the 
left-right action: GL('Hs) P = {g £ GL^Hs), gpg^ = p}- The group of phases by equivariance is U(Hs) x GL(T-Ls) P - 
But since we will use only phases by invariance for isospectral transformations (inner to (a) with <7 = 7 td(p)), 
the group is reduced to U(T-Ls )■ Because 3 g p £ GL(J~Ls) such that p = g P og\ (with a = 7 td(p), g P = s p u p 
with u p £ U(Hs) an unitary transformation and s p £ GL('Hs)/U('Hs) a SLOCC transformation), we have 
GLCHs) P = GL('Hs) g ag i = g p GL('Hs) a g~ 1 . The isotropy subgroups are then all isomorphic to the same model 
group for all p chosen in a single stratum (but it is different between two strata). 

To simplify the notations, we denote by K = U(7is) — U{n) (or U(Ha) if it acts on the right) the group of the 
phases by equivariance (or the group of ancilla transformations), and by G = GL{7is) — GL(n , C) the group of the 
non-unitary phases. For a stratum IT(%s) we denote by H l = GL{'Hs)(j — GL(qi) kl x ... x GL(qi) kl the stabilizer of 
the diagonal density operators with sorted eigenvalues (fc,; is the number of eigenvalues with degeneracy equal to qi). 
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The group of all phases by invariance of the stratum E l ('Hs) is the normal closure of H l : R l = (J geG gH l g 1 (H' is 
the smallest normal subgroup of G including H l ). 

We can also say that K is the group of the inobservable reconfigurations of the ancilla, G is the group of the unitary 
and SLOCC transformations of the system performed by the physicist, and H l is the subgroup of inefficient transfor¬ 
mations of the system (a transformation is inefficient with respect to a particular density operator). 

The non-unitary operator valued phases have been defined only with respect to the norm of the C*-module (the 
projection of purified states onto density operators). But the non-commutativity of the C*-algebra B{TLs) and the 
nonlinearity of the Hamiltonian of the purified dynamics (Hu{») = (H$ — ly^T^F*,) ® 1 Ha + 57*T fc 0 rj- (•)) induce 
a difficulty, because the phases do not commute with the generator of the dynamics. 

Definition 3 We call phase with respect to the generator of the dynamics, an operator g € G such that 

H u (g*)g* = gH u (*)V ( 81 ) 

For closed quantum systems where the Hamiltonian is linear and the phases are scalars, this condition is trivial. 

Proposition 1 The group of the phases with respect to the generator of the dynamics is Gc = Gfjeft H fj fe Gr k 
(H e ft = H$— | 7 fc rj.Tfcj, where the isotropy subgroups are defined for the adjoint action (G H eff ={j€ G, g~ 1 H e ^ g = 


H e ff). 

Proof: Let G B(TLs) be such that \F(IF^) = T. 

H u (g*)g* = H eff ® 1 nA g* + ^ k T k ® r* ($*)$* ( 82 ) 

= H eff g ® l Uj ySi + k T k g ® {{gW*YT k gW^ 4- (83) 

= H eff g ® 1«^ + l -l k T k g 0 ((Wj fl - 1 r fcfl W*) Tt 'F (84) 

H ef lg = gH eff if g G G H e ff , F k g = gT k if g G G Tk and g~ 1 F k g = F k if g G G Tk ■ H u {g^)g^> = gH u (^)^ if 
g G G H eff n fife Gr k - □ 


We remark that Vfc G K, 74/(l^ s ®fc T 'I , )l Ws ®A; T 4' = 1 Hs ®k T Hui'F)’^ (because r|(l^ s ®fc T 'F) = [k 1 WjfT k Wq,k) l: ^ = 

kwi^wr 1 )- 

We remark moreover that if the system Hamiltonian Hs and/or the jump operators Tfc are time-dependent, Gc is 
time-dependent. In that case, a possibility to avoid difficulties is to consider Gc = fl; Gx(t) but this group can be 
reduced to the unitary center of B(TLs) (U( 1) in finite dimension) except if Vi, Hs{t) and F k (t) belong to a same 
(small) subalgebra of B(TLs)- 


2. Operator valued geometric phases 

For closed quantum systems, the geometric phase concept is related to the cyclicity of the projected dynamics. Let 
[0, T] 3 i G p(t) G T>o(fHs) be a density operator solution of the Lindblad equation iHp = £(p). We say that the 
projected dynamics is cyclic if ttjj(p(T)) = 7Td(p(0)) (the spectrum of the statistical probabilities is the same at the 
start and at the end of the dynamics). We search a density operator [0, T] B t p(t) G VifHs) such that 

p(T) = p{Q) (85) 

\/t,3g(t)€G, p(t) = g{t)p(t)g(t) ] (86) 

p is the cyclic density operator associated with the density operator of cyclic projection. We note that g[t) G G and 
not K because there is no reason for which the transformation of p into a cyclic density operator implies no SLOCC 
operations. 

Theorem 3 A density operator p(t) with cyclic projection and its cyclic density operator pit) are related by p(t) = 
g(t)p(t)g(ty with the non-unitary operator valued phase g{t) defined by 

g(j\ — Tg - 1 ^ -1 fo J e ~ fo 


(87) 
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where the dynamical phase generator is defined by 

Vp € v(Hs), E{p) = Hs — ^ k r{T k + p k r kP rlp* e b(h s ) (88) 

and the geometric phase generator of the first kind is defined by 

% = dWW* ££1 1 {B(U s )M'Hs)) (89) 

iyy- th e inner product with W- viewed as a tangent vector, W p (t) being a standard purification of p(t): 1 = 

Wp(t)Wp(t)*. The geometric phase generator of the second kind is defined as 

77(1) = Wp(t)k(t)k(t)- l W~p(ty + h(t) e fj 1 (90) 

with k(t) G K and h(t) G B (ker p(t)) (t) 1 is the Lie algebm of H l with i the index of the stratum of p). 

We can note that 


2l||*(W0ll* = (*(W)\d.*(W)) t = tv nA \d*(W)))((*(W)\ (91) 

This definition is similar to the definition of the non-unitary geometric phase generator of non-hermitian closed 
quantum systems (with the C*-module in place of the Hilbert space of the system). It is moreover the non-adiabatic 
generalisation of the C*-geometric phase introduced in ref i 19 i 20 i 32 i 33 . We can also show that 

dp = + p^ (92) 

The definition of the dynamical phase E ( g P gX '> dt induces that the definition of total phase g(t) is implicit (the 

expression of g depends on itself). This is a consequence of the nonlinearity of the Schrodinger equation of the purified 
dynamics (we can say that the dynamical phase is a “nonlinear dynamical phase”). 

We can remark that h(t) = 0 if the cyclic dynamics ( 1 H > p(t)) takes place in regular strata. 

Proof: Since p = gpg ', 3k G K such that W p = gW p k. The nonlinear Schrodinger equation (for a standard purification 
state) is 

ihWp = H e ffw p + l - 1 k T k WpTlya(w p )y ( 93 ) 

= H e ffWp + l -T k T k Wp{W;V k W p y (94) 

= H eff Wp + l -l k T k W p WlY\{W;y (95) 

with H e ff = H S - ^ k T\T k . It follows that 

gW-pk + gWpk + gW-pk = -zh~ x (^H e ^gW p k + ^TtgWpWtg^ig^-^Wi^k) (96) 

and then 

g~ l gW p + Wpkk- 1 = -Wp - zh- 1 (g^H^gWp + (97) 

By multiplying on the right this expression by W~ we find 

g^g + W p kk~ 1 W~ = -W P W± - zh^g^Eigpg^g - h (98) 

with h G Bfkei p). We set 77 = W p kk~ 1 W ~ + h. Let g<% be such that g = Te~ lh E i g P gX ) dt g^_ 

9 = -ih- x E(gpgy g + Te "^" 1 £ E ^) dt g m (99) 

g~ x g = -ih^g^Eigpgyg + g^m (100) 

9^9* = -W p Wp — 9 => g<& = Te-fo(^w P +v)dt (101) 


It follows that 
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One needs now only to show that p G fj®: 



pp + pp' 

- Wpkk ~ 1 W\ + Wpk dk d W\ 

( 102 ) 


= Wpkk^wt - Wpkk^wt 

(103) 


= 0 

(104) 


since p = WpWj, = —fc 1 an d hp = 0 . Because of g~ 1 H l gp = GL(V.s)p = {g G GL(W. S ), gpg^ = p} 

(p = gpTTD(p)gl), the Lie algebra of GL(Hs)p is g[(' Hs)p = {X G B(Hs),Xp + pX t = 0}. We have then 

v e 0[(Ks)p c fjL □ 


From the point of view of the purified dynamics, we have 

T p (f) = Te - * R-1 Jo (i w^ a + ? 7 )dt 0 (i05) 


with T p = *(W P ), and « V p_2l||'Fp || 2 = (’Lpl'Lp)*- 

The right geometric phase kit) G K is arbitrary in the sense that nnsigWk) = gpg t = nnsigW), but it 

p 

induces the geometric phase generator of the second kind pk = Wpkk~ x W~ in the left geometric phase. Let 
A R G £l 1 (B(Hs),B('Hs)) be the generator of the right geometric phase: 

k(t) = Te~ f ° i *P«) ARdt ' (106) 

A possible choice consists to use the defintion of the Uhlmann geometric phas o 11 ' 12 : 

dW = uhl % W + WA R (107) 

where uhl % is solution of the following equation: 

dp = ™ 2 l p + p uhl 21 uhl< & ] = uhl % (108) 


This choice is set in order to W\\ = Wpk satisfies W|jW|| = W,jVF|| <*=> Wj|VF|* = which is the Uhlmann’s 

definition of the parallel transport for the density matrice s 11 ’ 12 . With this choice of right geometric phase, we have: 


21 + Wpkk- 1 Wi = i w _ (a + WA R W*) 

(109) 

= i w _ ( dWW * + WA R W*) 

( 110 ) 


( 111 ) 

The left generator appearing in the definition of the Uhlmann (right) geometric phase, 
geometric phase: 

generates then the total left 

'J'g - /o T 0Wp(t )^+ r lk(t))dt _ rjp e — / 0 T ®Wp(t) 1 hl< &dt 

( 112 ) 

= P e-fc™* 

(113) 


where C is the closed curve parametrized by [0,T] 3 t W p (t ) in B^Hs) — 9H riX ra( ( C) viewed as a manifold (Pe 

denoting the path-anti-ordered exponential). 

Since the total left geometric phase is automatically adaptated to the choice of the right geometric phase (by the 
presence of the left geometric phase generator of the second kind pk = Wpkk~ 1 W~), no natural choice of right 
geometric phase is imposed by the nonlinear Schrodinger equation and its associated Lindblad equation. Other 
definitions of the right geometric phase corresponding to other Uhlmann like connections can be choiced. This is the 
reason for which Dittman, Uhlmann and Rudolph^^ have find large classes of connections defining right geometric 
phases and compatible with the density matrix theory. 
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3. Special cases of left geometric phases 


a. Geometric phase with respect to the generator of the dynamics: We suppose that the dynamics takes place 
in regular strata. Let be the Lie algebra of Gc and let V'X' £ his ® "Hq, Sc( V I / ) = flcll/ be the Hilbert subspace 
which is the orbit of ’L by g^. Let {X,;}j be a set of generators of the Lie algebra g^ such that tr(X'dXj) = St. By 
construction {^(XjW)},; constitutes a basis of ^(^(W)) (with W £ B(hLs))- In spite of the orthonormalisation of 
the generators of g^, the basis {^(Xj W)}i is not orthonormal. Let {^(X^Wl)” 1 )}^ be the associated biorthonormal 
basis: 


«4'(X i (W t )- 1 )|^(X j W))) = (X i (Wt)- 1 |X j W) ffS 

(114) 

= tr (W^X^XjW) 

(115) 

= tr (WW^X^Xj) 

(116) 

= tr(X‘t Xj) 

(117) 


(118) 

(where we have used the cyclicity of the trace). Let ^f(Pc(W)) be the (non-orthogonal) projection onto <S(c(\lf(W)) 

in his ® hi A defined by 


PciW) = |X i W) ffS (X I (W t r 1 | = XiW tr (W-'X*.) 

(119) 

Let A L £ ^(BfHs), gc) be defined as 


A L \\V(W)\\t = (*(W)\*(P c (W))[d*(W)}}* 

(120) 

A L WW^ = P c (W)[dW}W f 

(121) 

A L WW f = tr (W^X^dW) X.;WW f 

(122) 

A l =tr(X if dWW~ 1 )X i 

(123) 

A L = tr (W ;t 2t) Xi 

(124) 


A L is the part of 21 which induces a phase with respect to the generator of the dynamics. 


Proposition 2 21 = i] y-A l (the geometric phase is a phase with respect to the generator of the dynamics) if and 

only if E(p) £ g £ ; 


[E(p),H e ff] 

= 0 


[E(P), IV 

1 = 0 


'e(p), if 

= 0 

W 0 k 

~ 1 — Wokk~ 1 . ’ 


Proof: i ^_21 = W P W P and W p = -g^gWp+g' 
i^„2l e gc if g £ G c and E(p) £ g £ . 

Note that g € Gc => E(gpg^) = gE^g -1 . It follows that 


(125) 

(126) 

(127) 

We have then *^,-.21 = ~g~ 1 g—ih~ 1 g~ 1 E(p)g—r]k- 


□ 


g{t) = Te~ lh 1 fo 9(t')E(p(t'))g(t') 1 dt'i> e ^o( l w. (tl) A +Vk)dt 
- ih _1 /* E(p(t'))dt'-f* iw ( A L + r] k )dt’ 

= Jie p( t > 


(128) 

(129) 


In that case, the nonlinearity of the dynamical phase does not then occur but the dynamical and the geometric phases 
are not seperated. 


b. Geometric phase by invariance: We suppose that the dynamics is inner to a regular stratum E l (Hs)- Let 
H l ~ GL(qi) x ... x GL(qi) be the stabilizer of the stratum E ’ l i(Hs) (we do not consider necessary that qj ^ qk for 
j ^ k). Let ff ~ u(gi) ® ... © u(qi) be its Lie algebra. Vcr £ E l ((Hs), let V a : B((Hs) —> ff be defined by 

V a {X) = P^XP* ( 130 ) 

where Pf is the orthogonal projection onto the ^-dimensional eigenspace of a associated with the eigenvalue p,,- (which 
has a degeneracy equal to qj). Let p £ V((Hs) be such that 7 td(p) = o. We know that 3g p £ G such that p = g P crgj, 
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and GL('H s ) p = g p H l g p 1 . Let V p : B(TLs) —> fll ('H-s)p be defined by V P (X) = P JP XP p with P p = g p P° g p 1 . 
Let A BL £ f)*) be defined as 

A bl = V ww \%) = pi wwX dWW~ 1 pW w ^ 

A BL is the part of 21 which induces a phase by invariance. 


(131) 


Proposition 3 iy^M l = i^ v _A BL (the geometric phase is a phase by invariance) if and only if E gp (a) £ ff with 
E gp {cr) = g~ 1 E(g p ag~ 1 )g p (g p £ U(U S ), p(t ) = g p (t)a{t)g p (t)- 1 £ D 0 (H S ))- 

Proof: i ^_21 = - g~ x g - ih~ 1 g- 1 E(p)g - g k and then i ^_21 £ gl(' H S ) P if g G GL(H S ) P and E{p) £ gl(" H s ) p <^=>- 
9p lE (p)d P e ()*• □ 

The geometric phase (and the dynamical phase) are phase by invariance, i.e. p{T) = p( 0) (the density operator is 
cyclic), if 


9p 1 H S g P ~ 7p k 9 P lr I r fcffp + 2 T9 P l ^kg P ^9p ^ k 9 P (r 1 g p 1 € ff 

The elements of gl {'Hs) P acting on the left can be converted as elements acting on the right: 

X£Ql(H s ) P ^Xp + pX t = 0 

=> XWW f + WW*X* = 0 

+WW t A t (W t )“ 1 = 0 
=> xw = —W{W~ 1 XW)^ 

The left action of X £ &l{'Hs)p is then equivalent to the right action of — {W~ x XW)^. We have then 


mV = (21 - A BL )W + WA 


\BR 


with 


A br = -(W 1 A SI 'W) t 

= -W* (p jWW 'y (W-^dWi (p^y (W _1 )t 


(132) 

(133) 

(134) 

(135) 

(136) 

(137) 

(138) 

(139) 


Pj vw f = gww^Pj 9wwi an( 4 W = gwwwf^i we have then W l P^ w] W = y/a 1 P° y/a = P!f since P° yfa = yfaP/ 
{P/ is an eigenprojection of a). We have then 


A br = -pi a dW ] (W ] )- l P° 


(140) 


Let W £ ShsB(Hs )■ dW^(W^) 1 = dyfayfa + y/adg^ w] g ww \yfa . Since da = dyfayfa + yfodyfa 


- = Oy/ay/C 1 1 1—1 1 1 

dyfayfa 1 = daa~ x — yfadyfaa^ 1 we have 

dWfWf~ x = daa~ l - yfadyfaa^ 1 - yfagf^ w ^dg WW ]yfa 1 (W £ ShsB/Hs)) 
But W^dWa^ 1 = yfagf^ w1 dgwwi 1 + yfadyfaa~ x , and then 

dW t (W t )- 1 = daa~ x - WUWa~ x (W £ S hs B{H s )) 

Finally, the restriction of A BR onto ShsB/Hs) is 

A\s R HS B(Hs) = P° a WUWP°a- x - P^daa~ x Pf 


(141) 


(142) 


(143) 


This is the expression of the generator of the Sjoqvist-Andersson geometric phased— (note that in the works of 
Sjoqvist and Andersson, da~ l = 0 since they consider only isospectral dynamics). A BR = —PPdWfW^)~ l Pf is 
then the non-unitary generalization of the generator of the Sjoqvist-Andersson geometric phase. We can note that 
equation (11371) can rewritten as 


dW = (21 - A bl )W + WA BR 


( 144 ) 
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which is an equation defining an Uhlmann like connection. In other words, we can choose the right geometric phase 
as being the Sjoqvist geometric phase in place of the Uhlmann geometric phase (fcfc -1 = — A BR ). In that case, the 
total left geometric phase is such that 

Te ~f?(iw ii m*+Vk(t))dt e g/g p (145) 

We will say that this Sjoqvist geometric phase is two-sided in the sense that it can be considered on the left ( A BL ) 
as a inefficient transformation of the system, or it can be considered on the right ( A BR ) as a reconfiguration of the 
ancilla. Its invariance nature implies that it can be revealed only by interferometry in accordance with its discovery 
by Sjoqvist etali&i£. 

c. Adiabatic geometric phase : Suppose that t H > Hs(t) and t H > Tk(t) are continuous and cyclic operators with 
respect to the time: H$(T) = Hs( 0) and Tfc(T) = 1^(0). Moreover suppose that |y| = max*, |-/ fe | <C 1. The nonlinear 
Schrodinger equation for the purified dynamics can be written as 

iM = H eff ® 1 H + ^7 fc r fc ® (IU (0)t r t fe IU (0)t *) T 4' + e>(| 7 | 2 ) (146) 

with 4M°) = be such that: 

= H s <8 1 u A y {0) (147) 

At the first order of perturbation, the nonlinearity is replaced by the pre-integration of the zero-order solution. 
Suppose that = 0) = ip(t = 0) 0 then 4/(°^(t) = t /A°)(t) ® with ihip = Hsip^ and 

W (0) = |^ ( 0 ) )(^a| and W (0) * = |^a)(V’ ( ° ) l (148) 


iu ( 0 ) *iw (0) = (V’ (0) |r fc |^ (0) )|? a )(Ca 


(149) 


Let A b ° (t) be the instantaneous eigenvalues of Hs(t ) (supposed non-degenerate) and C^\t) be the associated nor¬ 


malized eigenvectors: 

%<?> = 4V 

At the first order of the perturbation, the eigenvalues and the eigenvectors of H e ^{t) are 


Mb = m! 0) - ^7 fe <c£ 0 ) |r!r fc |c( 0) > + o(| 7 | 2 ) 


c b =cr^E' 


E 0) |rir fc |ci 0) ) 


2 U m1 0) - ^ - ^7‘((Ci 0) |r|r,|c[ 0) ) - (Cfirl^K™))-' 

The eigenvalues and the eigenvectors of Hu are then 

A b/3 = ^ ^7 fc <ci 0) I r fc I ci 0) ) (^ (0) |rt |^ (0) > + 0(l7l 2 ) 


Cl 0) +(9(|7| 2 ) 


(150) 

(151) 

(152) 

(153) 


^b/3 = C b ®£/J + k|J^7 


<CS 0) I r fc I c£ 0) ) <-0 (0) I rl-1 V’ (0) > 


d, 0) ®£ a + O(l7| 2 ) (154) 


- Mb - hid + |y((Ci 0) |rz|Ci 0) } - <C' 0) |r i |Ci 0) ))(V- (0) |rIlV’ (0) )- d 

where {£p)p is an arbitrary (time independent) basis of Ha- The eigenvectors of Hy are 




(Ci 0) irllCi 0) )^ (0) |r fc |^ 0) ) 


U ' Mb-M d -iy((Ci 0 ) |rjKi 0) }-(Cf|rt|cf})(^( o )|r^(°))^ )0ea + <!?(l71 } (155) 


with the eigenvectors of which are biorthonormal to the eigenvectors of H e ^ ((C*IC b ) = S c b). 


We suppose that 4/(t = 0) = $ aa (f = 0), A QQ being non-degenerate. If the evolution is slow and under some 
technical assumptions, we can prove the following adiabatic transport formula (see ref.-S) : 


pit) = g(t)p aa (t)g(tf + O ^rnax 


Yi\‘ 


(156) 
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g{t) = Te _lR_1 f ° E^{t)dt Te -1? A£\t)dt 


(157) 


with 


Paa = tr^l t aa »«t aa | = IC^XC^I + 0(171 


^ 1) = E A ^ic£ ) )(c£ ) * 


(158) 

(159) 


/)(!) _ Y^/a( 1 )*|_|a( 1 )\|a( 1 )\/a( 1 )* 
“ Vlba 1 dt Eca ' Eba 1 'hca 

be 


(160) 


c« = c + -V 

-Iba —b 0 Z-, 


<d 0) |r/c|d 0) X^ (0) |rllV’ (0) ) 


2 ^ ' Mb - w + |y«Ci 0) |ri|cD - (C' 0) |rdd 0) ))<V> (0) |r^( 0) )- d 


-( 0 ) 


(161) 


A(i)* = A _ l _ 
hba -lb 9 Z-, ' 


<cS° ) l r fclCf ) XV' (0) |r fc |^ (0) > 


2 ^ ^ > - (Q° l r /l^° ))<^ (0) l r # (0) ) 


( 0 ) 


(162) 


The adiabatic geometric phase issued from the adiabatic theorem^ Te f° A “ > ^ dt is related to the non-adiabatic 
geometric phase as follows. The role of p{t) in theorem [3] is played by p aa (t). We have W Paa = ICd)(£« I and 
W* = l^aXC^*!- It follows that 

fJaa. 1 ^ ' x —aa 1 


C 1 )\/a^C 1 )* I 


Hi, 2l = W p W* =1— c l ; )(C' 

W p aa Pact p aa l pf±aa , '2-i 


(163) 


and 


^ 1) = v».. a +Ei s Z ) >(Z ) 'i + o(iT) 

c^a 


(164) 


where we have use the fact that (C^iCi 1 ^*)b constitutes (at the order |y| 2 ) a biorthonormalized basis of Hs 
(J2b ICftaXiba*! = 1 )' More0ver 


+*»|ElIZ’XZ 1 '! ] = °(H 2 ) 


c^a 

d /dl)\/A(l)*| 




dt~ ca - ca 


(165) 


It follows that 77 = £ c#a e 0 l (H s ) Paa and then 


A « ] =i w 0 21 + V 


(166) 


The adiabatic geometric phase generator found in ref.— is the sum of the two geometric phase generators found 
theorem [3] with the “density eigenmatrix” p aa playing the role of /5(f). 

Let P, a = I—ha))«—ba I be the eigenprojection in Hs ® Ha associated with the eigenvectors of Hu related to £ a 

by perturbation. We consider the following reduced generator of geometric phases: 


( \r-AluJmj) A. = E(Z 1 ‘XlZ , >lZ 1 ><Z , 'l + °<M 2 > 


= l*00 + °<« , > 


(167) 


(168) 
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TABLE III. The different geometric phases of the open quantum systems. 


geometric phase 

generator 

values 

side 

type 

interpretation 

C *-geometric phase 

21 

G 

left 

by equivariance 

transformations to cyclic density operators 

Uhlmann geometric phase 

a r 

K 

right 

by invariance 

transition probability parallelism 

Sjoqvist geometric phase 

A bl / A br 

W 

two-sided 

by invariance 

interferometric phase 

reduced C *-geometric phase 

A L 

G c 

left 

w.r.t. the dynamics 

closest to the adiabatic phase 


*-\p. a )[w Paa }w; aa = 4 1 ) + 0 (| 7 | 2 ) (169) 

= i w Paoi % + V + 0( | 7 | 2 ) (170) 

Moreover 

E, i 1} ® |Q<£*| = P. a H u (^ aa )P. a (171) 

= E(p aa )®\Z a )(Z a \+0(\'y\ 2 ) (172) 

The geometric and dynamical phase generators appearing in the adiabatic theorem ref . 33 are clearly the adiabatic 
limit of the generators introduced in theorem [3] 

Let G , {$ 6 ( 3 } b p(i) C U(n s ) be the group of operators of Hs acting on the eigenvectors of Hu as phase changes. By 
definition we have 


e H u (g§> aa )g$ aa = A act5 l aQ + 0(| 7 | 2 ) (173) 

= 9Hu(^ aa m aa + 0( | 7 | 2 ) (174) 

Clearly E(gp aa g~ 1 ) = E(p aa ) and G^ bfj y bf) C Gc- The adiabatic geometric phase is a phase with respect to the gen¬ 
erator of the dynamics and A^ = A L . Finally the adiabatic assumption can be rewritten as 21 ~ iyy A L 
in accordance with the discussion found in ref.—. The geometric structure involved by the adiabatic geometric phases 
has been extensively studied in ref— and the interpretation of the adiabatic operator valued geometric phase has 
been studied in ref: 32 . 

This approach for the adiabatic geometric phases, which is usefull for bipartite quantum systems, is strongly limited 
for the open quantum systems described by a Lindblad equation. Indeed both the perturbative assumption and the 
adiabatic assumption are too drastic because they involve that the dynamics remains in the neighbourhood of the 
singular stratum of the pure states ( p aa is a pure state at the first order of perturbation). This excludes the more 
interesting cases where the relaxation rates are sufficiently large to deviate the dynamics or where the interaction 
duration is sufficiently large to the dynamics goes to the steady state (with the open quantum systems there is 
a competition between the adiabatic and the relaxation processes). In second order of perturbation the adiabatic 
approximation deals with mixed stated but the nonlinearity could induces strong difficulties to applied this approach. 
A more general approach of the adiabatic geometric phases could consists to use the notion of noncommutative eigen¬ 
values as introduced in ref— for bipartite quantum systems, but the nonlinearity could also induces strong difficulties. 

Table HI summarizes the different geometric phases of the open quantum systems. 


4. Example: control of a two-level system with dephasing and spontaneous emission 


We return to the illustrative example introduced in the previous section. Let t K > p(t ) be the evolution described 
by the left part of figure [3] By construction lim £ _).o ttd(p(0)) = linp_>. +(X) nD{p(t)) = S °(" Hs ) (the vertex associated 
with the pure states). The projected dynamics is then cyclic. The definition of the geometric phase depends on the 
choice of the local section t*->W(t)€ B(Hs) such that 3 g(t) € G, p{t) = g(t)W(t)W (t)t g(t)^ and lim t _> +00 W(t) = 
lining W’(O). We consider three interesting cases: 


the section y/a G T(£(' Hs),B{'Hs)) associated with the global section of the composite bundle: yja(t) = 


[ VW) 


0 


V 0 y/1 - p(t) 

section of diagonal density matrices); 


where t i-T p(t) = min Sp(p(f)) G 11(2) = [0, i] (a G T(T,(Hs),P{'Hs)) being the global 
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FIG. 6. Evolution of the populations and of the coherences of er(f) (left), p str ad (t ) (middle), and p wk ad (t ) (right). 



FIG. 7. Left: Representation of the dynamics t p(f), t i-> <r(f), 1 H> p atr ad (t) and 1 1 —> p wk ad (t) on TUHs) viewed as a ball 
(the spherical boundary is T>°{T-Ls) — CP 1 ~ S 2 the space of pure states, and the center is ^ 1 g 2 j/ 2 ) } the point of 

the microcanonical state). Each trajectory starts from the south pole, t cr{t), t i-» p str ~ ad (t) and t t-> p wk ~ ad (t) are closed 
whereas t e->- p(t) ends in the north pole. Right: same as the left figure but with a representation of r D('Hs) in accordance with 
its fibred structure, the height line of the cone is E("H.s) — [0, ~], the apex is the point of the microcanonical state D°°(Hs), 
the base is the space of pure states D 0 ('Hs) — CP 1 ~ S 2 (represented as a disk by azimuthal projection) and each circular 
sections are a fiber n] 3 1 (p) ~ CP 1 ~ S 2 . no : 'D('Hs) —> E(Ps) is then the orthogonal projection onto the height line of the 
cone. 


• the section yjp str ad £ associated with the strong adiabatic approximation^, i.e. p str ad = 

str — ad 

ti'w^|$ii})(($n| with $ n defined by eq. (11541) applied on the present system. We remark that £ , tr _ a a is a 
pure state. 

• the section yjp wk ~ ad g V('D('Hs),B('Hs)) associated with a weak adiabatic approximation where p wk ~ ad (t) = 

Pi (t)ti u A Hii (0» ((tit (t) | + P2 (t) | $ 21 (i))) (($2i (t) I ■ 

The sections are represented figure [ 6 ] and [7] 

We cannot use p str ~ ad since it is a pure state which does not belong to the same stratum that p (the strong adiabatic 
approximation is not satisfied for this dynamics). In the two other cases we have 

t(i) = '&(g E (t)g<n(t)W(t)k(t)) (175) 

with k(t) £ K the right geometric phase, g<&(t) = Te g G the left geometric phase and 3 s(t) = 

Te~ lh fo E (p( t )) dt the dynamical phase. The geometric phases are drawn figures [ 8 ] and [9] The effects of the 

A — 

geometric phases are represented figure [lU] 
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FIG. 8. Right geometric phases (top) and left geometric phase (bottom) with respect to the time, for the global section a(t) 
(left) and the weakly adiabatic dynamics p wk ^ ad {t) (right). We have use the decomposition of the operators onto {<ro, ax,a v ,a z } 
(ct 0 = lea) viewed as a basis of $Z 2 x 2 (C). 


IV. GEOMETRY OF MIXED STATES 

A. The geometry as a stratified principal composite bibundle 
1. Regular strata 

Let E*(%s) be a regular stratum, T> 1 (Hs) = (T, 1 (Hs)) be the stratum of the density operators over £*(%£), 

and B l (T-Ls) = n] 1 1 s ('Hs) be the stratum of the standard purified states over V l ('Hs)- 
Let L H g : G —> h\itB l (T-Ls) be the left action of G on B l ('Hs) defined by 

\/g£GyW £B l (Hs) L HS (g)W = gW (176) 

and Rhs ■ K — > Aut B t ('Hs) be the right action of K on B l (Hs) defined by 

Vfc e K,VW e B i {'Hs) R H s{k)W = Wk (177) 

These actions of G and K are free. The right action of K can be locally transformed as a left action of H l , indeed 

R H s(k)W = LnsiWkW-^W (178) 

Because iths(Rhs{ k)W) = ■KHs(LHs{WkW~ 1 )W) = tths(W) {{WkW~ 1 )WW^{WkW -1 )^ = WW^), we have 
WkW 1 6 GL('U s ) ww t C H i . 

B l (Rs) cannot be viewed as a left-bundle on 1 l’('Hs) since nnsigW) = g^HsiW)g\ but since ttd o nnsigW) = 
To ° khs(W) we can see B l (Hs) as a left principal G-bundle P' over Let <pp be the trivialization of PL 

4>p-.Gx E \U S ) ->■ B*(H S ) 

(g,a) ^ g^ 


(179) 
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FIG. 9. Same as figure [5] but with a space representation of the operators. 


P z is a trivial bundle because of the presence of the global section er i—^ fa. 

Since nHs{Wk) = nns(W), we can see B l ('Hs) as a righ principal AT-bundle Q l over V l (Hs)- Let cj )q be the 
trivialization of Q ! : 


x K -f B‘{U S ) 

(. P,k ) 


(180) 


Q l is a trivial bundle because of the presence of the global section p v /p. 

Let ttq : G —> G / H l be the canonical projection (n l G {g) = gH l ). Let (p' l D be the trivialization of V l (/Hs) viewed as a 
bundle of manifolds: 


: G/W x E \Us) -t V\U S ) 
{gH\a) gag f 


(181) 


Let 7: H l —>• G/H 1 be the projection defined by 7 Tg(h) = 7 t£j(/i) (H l C G) which is such that TT^(gH l g 3 ) = gH l . 
The complete geometric structure is then defined by the following commutative diagram: 


W x E \U S ) 

xid | 

G/H 1 x E ,; (' H s ) 


Gx E i {H s ) 

ttq xid 

G/H 1 x E l (' Hs) 
Pr 2 

s i ( Hs) 


B l (Hs) f——— VHUs) X A' 


™HS 

-> V^Us) 


^(Hs) 


Pri 

V l {H S ) 


Due to the presence of three floors ( B l (/Hs ), V l (/Hs) and TPi/Hs)) the bundle of purification is a composite bundle^ 36 . 
Due to the both left and right structures, the bundle of purification is a bibundle^I. We can then say that B l (TLs) is 
a principal composite bibundle. 
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FIG. 10. Evolution of the occupation probabilities of \E'(W / (t)fc(f)) (top) and of (g&(t)W (f)) (bottom) for W(t) = \Ja(t) 
(left) and W(t) = \Jp wk ~ ad {t) (right). 


The bundle can be rewritten with an arbitrary purification space Hs ® Ha with the left and right actions defined 
by 


\/g £ G,V'F £ H s <8> Ha, L*(g)i& = g® l Hj j& (182) 

Vfc e K,W £ H s ®H a , R^{k)^ = l Hs ® k^ (183) 

with L HS {g)^{W) = V{Ly{g)W) and R H s(k)^(W) = ^{R^,(k)^{W)). 


2. Singular strata 

Let Y± l (Hs) be a singular stratum charaterized by p\ = 0 with a degeneracy equal to n — to. A density operator 
p £ V l (Hs) can be defined by two kinds of data, PR anp the projection onto its support space and an order m invertible 
density matrix g £ 9Jl rnxrn (C) representing p in this space. Let (Xj)j=i,...,m be an orthonormal basis of Ranp. We 
have: 


P = 0 j k\Xj)(x k \ ^Ranp = \Xj){x J \ (184) 

with gi k = (x J | PXk)- Reciprocally we can introduce p : 2T(C m ) x V m (Hs) —*■ R l {Hs) defined by 

Vp e 2?‘(C m ),V( Xj )i= 1 ,...,m e V m (Hs), p(g, ( X j)j) = ^k\Xj){x k \ (185) 

where P ? (C m ) is the regular strata for the model Hilbert space C m such that E*(C m ) = 7 r£>( 2 ? l (C m )) is similar to 
T, z (Hs) deprived of the zeros. V m (H s ) = {(Xj)i=i,...,m S {Hs) m , (x 3 \Xk) = S{ } ~ R m (C n ) = u\n)/U(n-m) ( V m {C n ) 
is a Stiefel manifold). 

Let I a C {1, ..., n} be a set of indexes such that dist_ps(P a , PRanp) < f where P a = . gJa 10)(01 an( l dist_Fs(P, Q) = 
arccos | det(ZpZg )| 2 is the Fubini-Study distance on G m (Hs) = {P £ B{Hs),P 1 = P, Pi = P, dimRanP = to} ~ 
G m (C n ) (Zp is the matrix representing an arbitrary orthonormal basis of RanP in the basis (Cj)j)- There exists a 
basis (Xj)j °f Ranp such thafc2& 


3 Cj /- £ C, Xj — Cj ) ( CjkCk 
k(£I a 


(186) 
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U a = {P £ G m {T~Ls ) j dist Fs(P,P a ) < f } is an open chart of G m (/Hs) with the coordinates map 
Vm(Hs) is a non-trivial principal [/(m)-bundle over Gm{T~Ls) with local trivializations 


(^yn{n—m) 

i c jk)jel a ,k£I‘ 


r s :U a x U(m) -»■ V m (H s )\u° 

( p ,u) * (z a P u)\ |0)<x fc l 


(187) 


with Zp £dJl mxm (C) be such that (Zf,)i k = (OHx'fc) (C Xj) jzia is an orthonormal basis of RanP). 

Let H i0 and H il be such that H l = H i0 x H a with H l ° ~ GL{m,C)a (a £ £*(C m )) and H u ~ GL(n - 
to,C). G/7P 1 ~ y m (C") = GL(n,C)/GL(n - m,C) = V m (C n ) X;y( m ) GL(m,C ) = {[(xw _1 ,uff),u £ f/(m)],x £ 
V m { < C n ),g £ GL(m, C)} is a non-compact Stiefel manifold. Let <jj l g be the local trivialization of G/H 1 1 induced by 
$j° : t/° x GL(m, C) -f V ro (C n ) 

(-P,5) [(Z^u _1 ,u5),u£ C7(m)] 

The complete geometric structure is defined by the following commutative diagram: 


H l x S \H S ) 
G/H i x H i0 x S '(Us) 

Pri X Pr 3 | 

G/H i x T.y-Hs) 


{0g a }aXid 


-4 GxS'(Us) 

7Tq xid I 

-4 G/H a x S i (H s ) 


-4 BtiUs) — T> i (C rn ) x Vm(Hs) 4 


idx{</>|}„ 


G/H' x STUs) 


"HS J 

-4 T> i ('Hs) 4- 


idXTTg J 

V^C” 1 ) x G m (Hs) 


V i {C m ) x Gm(H s ) x H(m) 

Pri X Pr 2 

X>'(C m ) x G m (Hs) 


Pr 2 j 

E*(«s) 


STHs) 4- 


7T D oPri J 

E'(C m ) 


where 7 t£j : G —> G/H 11 and ttq : G/H l1 —»• G/(H l1 x H l0 ) are the canonical projections; p*(g,P) = p(g,(Xi)i) f° r 
an arbitrary orthonormal basis (xi)i of RanP and 


P* ■ ^(Tis) x V m (Hs) -»• ^(Ws) 

(ftx) ^ v^lxWl 1 ; 

B'CHs) can be viewed as a (non-trivial) right principal [/(m)-bundle Q z over D‘(/Hs ) with local trivializations p* o 
</>S op/ 1 . Note that [/(to) acts on B l (/Hs) by the right action: 

Vu £ U(m),VW = p*(gwA a (gw,uw)) £ B l (Hs), Ri oc (u)W = p*(gw,(t>s(gw,uwu )) (189) 

By construction Vu £ U(m ), 3fc Ui w £ K such that Ri oc {u) = RHs(ku,w)- We denote by K'{, v the subgroup of K such 
that Rioc(U(m))\w = Rhs(Kw). On the left, B l (/Hs) is a trivial bundle P* over Y, l [/Hs) with typical fiber G/H ll 
(which is not a group since TP 1 is not normal). 


B. The geometry as categorical principal bundles 

We want to endow the geometric structure with a connective structure describing the geometric phases. On the 
regular strata, the geometric structure is a principal composite bibundle. A connective structure on a principal 
bibundlc 3 ^ and a connective structure on a principal composite bundle^ present a higher degree of complication than 
for an usual principal bundle, the compatibility between the two approaches can be also a difficulty. In order to avoid 
this, and to unify the geometric description between the left and the right, we propose to use a categorical geometric 
generalization of the principal bundle^ -25 ’ 39 . 


1. Left principal categorical bundle 

We consider a stratum Yi l (/Hs)- Let & 1 be the (trivial) category with the set of objects Obj^ 1 = E l CH 5 ) and the 
set of morphisms Morphs* = {id CT , a £ Let SZ l be the category with Obj 2?’ 1 = B l (/Hs) and Morphs = 
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{(W, k),W € £ K{y} (K’ w = K on the regular strata) with the source, the target and the identity maps 

defined by 


V(W, k) £ Morphs®, s{W, k) = W 

(190) 

t(W, k) = R H s(k)W = Wk 

(191) 

idw = ( W., 1 u s ) 

(192) 


the arrow composition being 

MW £ B l (H s ),Mk £ K z w ,Mk' £ K l wk , ( Wk , k') o (W, k ) = (W, /cfc') (193) 


By construction, arrows of different strata are not composable. Let wd £ Funct(^®, S fil ) be the full functor defined 

by 


MW £ Obj^‘, wd(W) = ttd o tths(W) (194) 

M{W,k) £ Morph^, w D (W,k) = id n D o^ HS (w) (195) 

Note that WD(t(W, k)) = wuiWk) = nD{Wkk~ 1 W 1) = t:d(WW^) = t(wn(W, k)). 

We can remark that the arrows could also be defined by an action on the left: Morphs* = {(h, W), W £ B l (kLs),h £ 
WK\yW* C H *} with t(h, W) = Lns(h)W, since 3k £ I< such that h = WkW* and Lns(WkW*)W = WkWW* = 
Wk{ 1 ^ P ker W ) = Wk (fcPker W = 0, Mk £ K l w ). _ 

Let M? 1 be the groupoid defined by ObjSf* = G and Morph^ 1 = G x iP with the source, the target and the identity 
maps defined by 


M{g, h) £ MorphSP, s{g , h) = g 
t(g, h) = gh 
id s = (g,l-H s ) 

the usual arrow composition (called the vertical composition of arrows) being 

Mg £ G, Mh, h' £ iP, (gh, h') o (g, h) = (g, hh') 

and the law of the semi-direct product of groups (called the horizontal composition of arrows) being 

M(g, h), ( gh') £ G x H\ (g', h')(g , h) = (g'g, g~ 1 tigh) 

Let 2zf : MorphSP —> EndFunct(^? 2 *) be the left action of ( S l onto M? 1 defined by 

M(g, h) £ Morph 9*,MW £ Obj^\ JV g , h (W) = L HS (gh)W 
M(W,k) £ Morphs, J^ g ^(W,k) = (ghW,k) 


or similarly by 


W 

( W,k ) 

Wk 


ghW 

( ghW,k) 

ghWk 


The composition of left actions is covariant with the horizontal composition of the groupoid arrows: 


■^g' ,h 




g,h 


= <£ , _ 

9 9,9 


1 h' gh 


— -^(g’ ,h’)(g,h) 


(196) 

(197) 

(198) 


(199) 


( 200 ) 


( 201 ) 

( 202 ) 


(203) 


( 204 ) 
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The compatibility between the projection functor vjd and the left action endofunctor Jz? is shown by the following 
commutative diagram (where aw = t^d ° ^hs{W)) : 


ghW 



(205) 


a w 

is naturally equivalent to Sf* x by the trivialization functors (j) 1 ^ G Funct(Sf* x y i , S? 1 ) and G 
Funct(,^*,Sf* x y i ) defined by 


(5,0-) 

(g,h, trade) 

(gh,a) 

(w yjow* ,a w ) 
{W^aw* ,h,&w ,i A„ w ) 

Wkyjaw* 


(gVo,k) 

g^/ak 

W 

(W, k) 

Wk 


with h = s faksfa' and aw = To ° ^HsiW)- We can note that (f> l ^,(g, a) = <f>p(g, a). 


(206) 


Viewed as a left principal categorial bundle, the presence of the both left and right actions takes a natural meaning. 
The right action of K (inobservable reconfiguration of the ancilla) defines the arrows of the purification category 
whereas the left action of G (unitary and SLOCC transformations of the system) defines endofunctors of the purifi¬ 
cation category playing the role of gauge changes. 


2. Right principal categorical bundle 

It is possible to reverse the description, by considering that the right action of K defines endofunctors of gauge 
changes (in accordance with their inobservable character) and that the left action of G defines the arrows. 


Let ./#* be the category defined by Obj.-#* = V l {Hs) and Morph./#* = Gx V l (Hs) 1 with the source, the target, 
and the identity maps defined by 


\/{g,p) G Morphs#*. 

. s{g,p) = p 

(207) 


Oi 

O. 

II 

3 

3 

(208) 


idp = (1 HsiP) 

(209) 

the arrow composition being defined by 



VpG2>*(Hs),V 5 , 5 'GG, 

(s', W') ° (g,p) = (, g'g,p ) 

(210) 

Let 3 l be the category defined by ObjcS* = B l (TLs) and 
the identity maps defined by 

MorphcS* =GxB i (TLs), with the 

source, the target and 

V(</, W) G MorphiF, 

s(g,W) = W 

(211) 


t(g,W)=gW 

(212) 


id w = (1 w s , W) 

(213) 
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the arrow composition being defined by 

WW€B i (n s )yg,g' €G, ( g',gW)o( g ,W) = (g'g,W ) (214) 

Let whs G Funct(=S*, be the full functor defined by 

VlFGObj=2 4 , whs(W) = t:hs{W) = WW ] (215) 

V(< 7 , W) G Morph=S*, wns{g, W) = (g, WW^) (216) 

Let JT be the groupoid with ObjJT = A"' and MorphJT = AT i x K‘ (where K\ v ~ A'', VIF G B l (H s ); K l = K 
on the regular strata), with the source, the target and the identity maps defined by 

V(g, k) G MorphJF 4 , s(q,k)=k (217) 

t(q,k)=qk (218) 

idfc = (lw s ,*0 (219) 

the usual arrow composition (called the vertical composition of arrows) being 

Vfc G AT*,V< 7 , </ G AT®, (, q , ,qk)o(q,k) = (q'q,k ) (220) 

and the law of the semi-direct product of groups (called the horizontal composition of arrows) being 

V(<7, k), (<?', fc') G iT K K\ (q’, k’)(q, k ) = {q’k’q{k’)-\k'k) (221) 

Let £% : Morph.yF 1 —> EndFunct(=2 l ) be the right action of JF 4 onto =2* defined by 

IF IF(/fc 

(ff.VUqfc) (222) 

gW gWqk 


where to simplify the notations, for the singular strata, we write Wk in place of Ri oc (k)W. The composition of right 
actions is contravariant with the horizontal composition of the groupoid arrows: 


&q,k ° &q',k' &q'k'q(k')- 1 ,k'k ^(q',k')(q,k) ( 223 ) 

The compatibility between the projection functor whs and the right action endofunctor £% is shown by the following 
commutative diagram (where pw = WW'): 


Wqk 


IF 


gWqk 


t g,Wqk ^ 


(9 


•H. 

gW 


Pw 

(a,pw) 

gpwg ] 


( 224 ) 


pw 

(9,Pw- 

gpwg 1 


i? 4 is naturally equivalent to ..<#* x by the trivialization functors (j) 1 ^ G Funct(./#* x JF*, J2 l ) and <} G 
Funct(=S*, x Jff 1 ) defined by 


(P, k) 

(g,p,q\fc)| 

(. gpg\qk) 

(1 wwWww t* IF) 

(g, WW t, yjWV* g y/wwf , vdFWT* W) 


sfpk 


apa^a\fp*^\fp k ) 




( ff TFlFt ff t ;V ^wWV 5 if) 


Vgpghk 

IF 

(g,w) 

gW 


( 225 ) 
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TABLE IV. The different entities in the category formalism. 


entities 

symb. 

left role 

right role 

interpretation 

left transformations 

9 

object gauge changes 

arrows of the categories 

unitary/SLOCC transformations of S 

right transformations 

k 

arrows of the total category 

object gauge changes 

inobservable reconfigurations of A 

stabilizers 

h 

arrow gauge changes 

inner arrows 

inefficient transformations of S 

diagonal matrices 

G 

objects of the base category 


statistical probabilities 

density operators 

p 


objects of the base category 

quantum mixed states 

bounded operators 

w 

objects of the total category 

objects of the total category 

quantum purified states 


TABLE V. The different category structures. 


Category 

symb. 

objects 

arrows 

interpretation 

left base category 

yi 

£*(' Hs) 


space of the statistical probabilities 

right base category 

JG 

V^Hs) 

G x V\U S ) 

quantum mixed state space with system transformations 

left total category 



B'(Us) 

quantum purified state space with ancilla transformations 

right total category 

& 

B*(' Hs ) 

G x 

quantum purified state space with system transformations 

left groupoid 

yi 

G 

GxIF 

unitary and SLOCC transformations of the system 

right groupoid 

jr 

K l 

K i x K j 

inobservable reconfigurations of the ancilla 


We can note that 4> l ^{p, k) = </>q (p, k ). The arrows of the groupoid are K l k K l because (\/ gpg 1 g^/p, k ) models the 
transition from to A' 2 ^ fc over the arrow ( g , p) £ Morph.-# 2 . 

We remark that we can unify the categories: & = |_ | £ and = |_|, Q l , on which act the groupoids = | _ | 2 Sf 2 

and Crff = |J J^ 2 (except that the arrows of different strata cannot be composed neither horizontally nor vertically). 
& and i? can be viewed as kinds of categorical principal bundles, but not strictly as principal 2-bundles2i^^ (the 
direct categorical generalization of the principal bundles) since in contrast with Sf 2 and J^ 2 , Sf and J# have not the 
structure of Lie crossed modules ^ 25 : 39 . 

Table IIVI summarizes the roles of the different entities in the category formalism, and table 0 summarizes the 
category structures. 


C. Connections 

1. The connective structure of the left principal categorical bundle 

A connection on a principal 2-bundle with base space being a trivial category is^Ir— the data of a familly of 
connections on the objects related by kinds of gauge change associated with the arrows. 

More precisely, let B l (fHs) = Obj^ 12 be viewed as a principal bundle over E 2 (77s). VW £ B l (jHs), the (left) vertical 
tangent space at W is defined by 

V£B i (H s ) = {XW,X e Q \t) a } (226) 

where g is the Lie algebra associated with G and t) 21 is the Lie algebra associated with H zl (f)* 1 = {0} on the regular 
strata). We define the (left) horizontal tangent space at W by 

H^B\Us) = {%£ T w B i {Hs), XW* = 0} (227) 

On the regular strata, H^B l (fHs) is a connection for the G-principal bundle B l (fHs)- This connection is characterized 
by the connection 1-form £ fl 1 (B l ('Hs), fl) defined by = XW~ x . On the singular strata, we consider the 

trivial G-bundle G x E 1 (Hs) and we set an arbitrary connection H[ g ^(G x T, l ('Hs)) such that x 

= Hk a )B x ('Hs)- We introduce £ kl 1 (B l (Hs), q) such that TTQ*(j>pUJ ^ is the connection 1-form of 

G x Two different choices of arbitrary connections in CT )(G x Y, l (Hs)) are related by = gw £ 

fl 1 (B l (Hs ), l)* 1 ) which is a fj 2l -gauge change of the second kind. 
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Let $ : Obj'P* —» G be a left equivariant map from to G (VIU £ Vg £ G, $(gIU) = g$(W)). $ 

defines a gauge transformation a$ : Obj^ 14 —>■ Obj,^ 4 by a$(IU) = Under this gauge transformation, the 

connection 1-form becomes: 

a $ w vv(-^) = (228) 

= c 4 {w)w (<S>(W)X) + to^ w)w (d^(WMW)- 1 ^(W)W) (229) 

since Vy(t) a curve on B l (%s) having X as tangent vector at t = 0 ( 7 ( 0 ) = W), we have Jj-a$( 7 (t)) = ^j($( 7 (i)) 7 (i)) = 
ft) + ^ftlftftu^’ and then a $*-£ = ^a*(7(i))| t=0 = d<f>$ -1< f>TU + (d<I><I> _1 : B' l {%s) -» 0 ). It follows 

that 

a%u^{X) = ^(W-Ou^X^IU)- 1 + ^(W^W)- 1 (230) 

Let T : Obj,^ 4 —> UmeB'(« s )-^vv be a right equivariant map (VW £ B l (Hs), V/c £ Y(Wfc) = T(VU)£;). 
T defines another kind of gauge transformation ar ■ Obj'P 4 —>• ObjU 4 by ax(IU) = WT(W). Under this gauge 
transformation of the second kind, the connection 1-form becomes 

ClyOJ^(X) = &wy(w) ( a Y*X) (231) 

= “wr ( w)(X T(W0) + ^ r{w) {WT{W)T{W)~ 1 dT{W)) (232) 

since £ar( 7 ft)) = ^^T(q(t)) + 7 (t) dT( 'J t (t)) and then ar*X = XT + WT~ 1 dT (T~ 1 dT : B’ l {Hs) ->■ 8^5 being 
the Lie algebra of K^). It follows that 

axu4(£) = uj^(X) + WdT{W)T(W)~ 1 W* (233) 

(IUdTT- 1 IU*)IUIU t + IUIU t (IUdTT- 1 IU*)t = WdTT^W+WTdT^W = 0, we have thenJUdT(IU)T(IU)^ 1 IU* £ 
0l(Hs)wt C f) 4 = © g eG The gauge transformation of the second kind is then a f) 4 gauge transformation. 

The connection on S? 1 is then not consituted by a single usual connection of B l (Hs ), but by a familly of usual 
connections of B l {7is) related by additions of (f-valued 1-forms. 

We are now able to provide the local expressions of the connection. Let s £ T(T, l ('hs),B l ('Hs)) be a local section 
(s:<74 W s (a)). We have 

s*lo l = dW s W* = s 2l £ ^(2*0 H S ),S) (234) 

The G-gauge potential s 2t of the connection is the generator of the operator valued geometric phase (the connection 
having being defined for that purpose). Under a gauge change of the first kind W s (cr) —> g(a)W s (a) (with g £ 
n°(S 4 CH 5 ),G)), the gauge potential becomes: 

° s % = g s %g- 1 +dgg- 1 (235) 

Under a gauge change of the second kind W s (cr) —> W s (a)k(a) (with k £ f2 4 (E 4 (%s), J), the gauge potential 

becomes: 

sfe 2l = s 2( + s ?7 fc (236) 

with s r]k = W s dkk~ l W* £ ft 1 (T, 1 (Hs), f) 4 ) (the f) 4 -potential-transformation becomes under a gauge change of the first 
kind gs r]k = g s Vkg~ 1 )- The arbitrary nature of the f) 4 -potential-transformation is clear in the formalism (it is gauge 
change of the second kind). 

The composite nature of the bundle B\T~Ls ) permits also to consider an intermediate entity between the local and 

t t 

the global expression of the connection. Let s* £ T (T> 1 {Hs) , B l (Hs)) be the local section (s*(p) = W t (p)) such that 

S0 

t 

S*(p) = ^hs(s(ttd{p)))- The 1-form: 

S 1 *UJ L = dW f W* = si 2l £ QftDft' Hs),g) (237) 

si 

satishes for all h £ Vt° (V 1 (T~L S ), H l ) such that h{p) £ GL(Hs) P = g p H l g ~ 1 {4>n(g P ,^D(p)) = P ), the following gauge 
transformation rule: 

hsi Ql = h si Qlh-l + dhh ,- 1 


( 238 ) 
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Property 3 \/g £ we have 

= 3 s *21(p)3 _1 + dgg _1 + g s p(g, p)g~ l (239) 

where s p £ f2 1 (G' x T> l (fHs), b l ) is a gauge change of the second kind. 

t 

Proof: Vp £ Hs ), 3/c p £ K such that s*(p) = W T = y/pk p . In a first time, we consider the case of the global 

t 

section s*(p) = yfp. 

dyfpyfp + yfpdyfp = dp =t> d^fp^fp + (d v / p v / p) t = i(dp + dp t ) (240) 

=> 3X p , s.t. Xj, = -X p , dy/p^/p=^dp + X p (241) 

dyfpy/l? = \dpp* + X pP * (242) 


t I" 

We denote by = d^fp^fp* the gauge potential for the global section, and by ld *2l(p) = \dpp*. By noting that 

, T . ,t —- 

X p p*p + p{Xpp*y = 0, we see that v *21 and ld *2l are related by a ff-gauge change. 
v /7, 2l(ppp t ) = ld *2l(ppp t ) + X gpgi g^ 1 p*g~ 1 

= l;d{gpg ] )(gpg ] Y + ^ gpg tP t_1 p*p Jl 

= \gdpp*g ~ 1 + dgg^ 1 + g Qpdptp t_1 p* - Tjd^dg + ff _ 1 X gpgt p t_ V*^ ST 1 
= g'^*QL{p)g~ 1 + dgg- 1 + g QpdpV~V - + S _1 ^ gPg tP t_1 p* - -YpP*^) S _1 


m(g,p) 


pip + pp I = 0, pi is a h*-gauge change. 
We return to the general case: 


s *2l(p) = dW fW* r = v^a (p) + Vpdkpk; 1 ^* 

s • s# 

*2t(5Pfl t ) = ff vC *2l(p)5 _1 + rfflP" 1 + PPi(s,p)p _1 + V QPg^dkgpgt k~ pgi \/pppt* 


(243) 

(244) 

(245) 

(246) 


(247) 


(248) 


= g sm< &(p)g 1 + dgg 1 + gpi(g, p)g 1 + \/gpg^dk gpg] k g ^ \Jgpg t - gjpdkpkp 1 yfp g 1 (249) 


9 V 2 (g,p)g 1 


P2P+Pp\ = g 1 \Jgpg^dkgpgi k gpg t \Jgpg t gp—y/pdkpkp 1 Vp+pp 1 ' \Jgpg t fc gpg t dk gpgi -Jgpg^g ] 1 -^pk p dk p 1 ^fp (250) 

Since kdk _1 = — dfc/c -1 and pg^ = g- 1 \Jgpg i \Jgpg i => pg^\fgpg' = p -1 a/ gpg\ we have P 2 P + pp\ = 0, 772 is a 
ff-gauge change. □ 

s *2l plays the roles of a H % - gauge potential and of a G/lP-connection 1-form. The local and the intermediate 
expressions of the connection are related by 

s *2l(p) = gp S $l{TT D (p))gy 1 + dg p g~ p 1 + g p ^p(g p , a)g~ 1 (251) 

Conversely we have 

s 2l(a) = ( 7 T* D si %)(a) (252) 

s 2l defines the curvin g 19 : 21- — of as 

s B l = cf 21 - s 2l A s 2l € Cl 2 ('E i ('Hs),W : ) 


( 253 ) 
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Property 4 s B L (cr) G ff 


Proof: Since SIR 7 10 = dWW t we have d2lWT0 — SI A d(WW^) = —dW A gH 0 . It follows that (MIUI0 — SI A 
StR 7 R’t - SI A = _ dW A d^t. We have then B L WWt = SI A WWW + dW A d_I0 and 0105^ = 

-2U00 A Sit + dW A dW*, implying that B L WW*< + WW f B L ^ = 0 => B L G 01(70) C 0 □ 

Remark, for a regular stratum, P 1 is flat: S I3 L = 0 ( s Sl = dW s W~ 1 is pure gauge). Under a gauge change of the 
first kind, the curving becomes 


gs B L = g s B L g ~ 1 


(254) 


and under a gauge change of the second kind, it becomes 

sk B L = s B l + d s m - “rjk A s r? fc - [ s 3l, s r, k } (255) 

As in refJ^^— we can define fake curvatures of tf* 1 , as being S F L = d s A L — S A L A S A L — S B L G fI 2 (E*(?0),0c) 
( S A L = tr(A0 s 3l)Xj, (Xi)i being generators of q c ) and S F BL = d s A BL - S A BL A S A BL - S B L G ti 2 {Y?(U s ),V) 
( S A BL = •pp( s 2l), P p being the projection of 13(70) onto fll(70)0 S F L is the fake curvature of PA 1 with respect 
to the generator of the dynamics and S F BL is the fake curvature of SZ 1 with respect to the Sjoqvist-Andersson 
connection. For the adiabatic case, interpretations of B L and F L can be found in refill, the curving is a measure 
of the “kinematic decoherence” (decoherence induced by variations in £*(70) during the dynamics), and the fake 
curvature is a measure of the non-adiabaticity of the system entangled with the ancilla. To understand these quantities 
in a non-adiabatic context, let ’F G r(£*(70),70 <8> ’Ha) be a local section of a singular stratum i. We denote by 
p = (0,..., 0,pi, G £(ra) (m < n = dim 70, p,; G]0, 1[) the local coordinates of cr G £*(70). By a Schmidt 

decomposition^ we can write that 


30(p) G 70,3** (p) g H a , ^{p) = '52y/Pi<t ,% (p)®X t (p) (256) 

2=1 

After some algebras, we can show that the curving associated with this section satisfies 


tr ■H s (p(p)B L {p)) = Y,MFa(p) ~ A s (p) A A s (p))u (257) 

2=1 

with p(p) = 7r^(^(p)); Fa{p) = dA^(p) + A^(p) A A^(p) G fl 2 (£*(70), 9Jl mX m(C)) is the curvature associated 
with the Berry potential of the ancilla vectors of the Schmidt decomposition A^(p)ij = (x l (p)\dx : ’ (p))ha- Except 
non-abelian corrections, the statistical average of the curving lo p (B l ) is essentially the average of the ancilla Berry 
curvature for the Schmidt decomposition. Let a a ,a b ,cr c G £*(70) be three infinitely close points. Let ( a a a b a c ) be 
the infinitesimal triangular simplex defined by these points. We have 


e -//(.“.F a _ z bc z ca z ab (258) 

with (. Z ab ) i:j = (Xi(p a )\Xj(p b ))• Now, if {x*(p 6 )}i= spans the same subspace that {x*(p a )}*=i,...,m> then Z ab is 
just a matrix of basis change inner to this subspace and | det(Z ab )\ = 1. If the subspace spanned by {Xi(p b )}i=i,...,m 
is different from the subspace spanned by {xi{p a )}i=i,...,m, then | det(Z ab )\ < 1. If follows that det(e _ df^a^b aC) Ua^ 
is a measure of the change of subspace in the neighbourhood of the simplex ( a a a b a c ). It follows that Fy\{p) measures 
the propensity of the dynamics to leave the subspace spanned by {Xi{p)}i=i,...,m- Finally the curving B L (p) measures 
the propensity of the dynamics to leave the subspace spanned by {(j>i(p) <S> Xi(p)}i=i,...,m- In order to describe the 
dynamics in the neighbourhood of p (denoted V(p)), we need a subspace spanned by Up'<=v(p){0(7 ,/ ) ® Xi(p')}i-i, 

But even if V(p) is infinitesimal, if ||I3 i (p)|| is large, the dimension of this subspace will be larger than m. Finally 
we can understand B L (p) as the measure of the propensity of the system to leave its initial singular stratum to a less 
singular stratum. In accordance to this interpretation, B L is zero in regular strata. With similar arguments F L and 
F bl measure the propensity of the system to leave Sc (\F) and to leave a regime involving only phases by invariance. 
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2. The connective structure of the right principal categorical bundle 


Let the following commutative diagramms: 


Obj& 


id 


■* Morphi?* 


Obji2* 4- 


Obj .#* 


id 


4 Morph.# 1 


WHS 

Obj.#* 4- 


Morphi?* 

ttHS 

Morph.#* 


Objig* 4- 


Morphi2* 


Obj.#' 4- 


Morph.#* 


where we consider 0bji2* as a /^-principal bundle on Obj.#* and Morph=3* as K 1 kK 1 -principal bundle on Morph,#*. 
A connection on a principal 2-bundle with a base space being a non-trivial category is^ the data of connections 
HwObjJS 1 and H^ g< w) Morphi?* compatible in the sense where 


id*-Hvy0bjc2* = H( lu ^ W )M orph=2* 
t*H( g ,W) MorphJ?* = IL g iy0bjc3* 
Morph=S* = iijyObji?* 


(259) 

(260) 
(261) 


the vertical tangent spaces being defined by VwObj£! 1 = V^B l {Hs) = {WX, X £ t l w } and Vjg^MorptuS* = 
{(0 ,WX),X £ £\y}- The connections are characterized by connection 1-forms, u> R £ 0 1 (0bj.#*,f*) and £ 
fi 1 (Morph.#*, t* 5 t*) related by 


id*w* = u R 

= V(tO R ) 


(262) 

(263) 

(264) 


7r { : r 5 e 


J &2 yXie . fci ^ hi _v hi 

and * ‘ . Let <; 0 £ F(Obj.#*, Obji?*) and ?_>. £ lYMorph.#*, Morphi?*) 

Y 3 X A Y 3 X Y + X v J ; ; 


with 

be the trivializing local sections: 


These local sections satisfy: 


So(p) = WpVh s ) = Vp 
s^{g,p) = <t>%(g,p, i u A Vu A ) = (Vgpg^Vp* >Vp) 

id, o{p ) = ?->(id p ) = (1 h a ,Vp) 

s{g^{g,p)) = So(s(g,p)) = Vp 


t(s^(h,p)) = s 0 (t{g,p)) = Vgpg ] 

The local data of the connection, the gauge potentials A R = s*w R £ f] 1 (Obj#*, 6*) and 
AV = £ fl 1 (Morph#*, t* 5 {*) are such that 

id* A* (p) = A R (p) 

^{A R {g,p)) = A*(p) 
tV e (AK(g,p)) = AK(gpV) 

These conditions imply that A R ^(g,p) = A R (p) + A R .(g,p) with A^(g,p) £ fi 1 (Morph#*,£*) such that 

A*(l„ s ,p) = 0 

A*(g,p) = A^(gpg^) - A R (p) 


(265) 

(266) 

(267) 

(268) 
(269) 


(270) 

(271) 

(272) 

(273) 

(274) 


In order to the horizontal lift of the arrows corresponding to the parallel transport of p induces the geometric 
phases, we must have 


(Pi l u A ) 

(g^pgliAnV 


W) g% \jg%pglV R 

I (sa,ffa 

(g*pgh,qk) \fmp^k R 



X(q,k) 

, 1 


, 

- > 

(g< 2 i,p,q,k) 

1 
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T 

_ r* i . Uhu^^/ _ rt ■ A R dt' 

where g% = Te Jo vp is the left geometric phase, k K = Te Jo Vp is the (right) Uhlmann geometric phase, 

and (q, k) G MorphJU* is the searched horizontal lift: 


(275) 

(276) 


/ 7 \ nn — fn i- j-A R dt' 

(q,k) = Te Jo —* 

/rjn — --A^dt' rm — f* 1 .-A^dt's 

= (Te Jo 9 2i >'/p 5 Te JQ vp ° j 

where i^A R = (Te ,vp A ^ dt i^-A Rr fe ^ 0 l B<n,s/p A ^ dt by virtue of the intermediate representation theorem^. 

We have in the definition of the parallel transport \Jg^pg^k R and not g<&^fpk R as for the right connection, because 
we want that the target of parallel transported arrow correspond to an usual parallel transport in the “target bundle”, 

i-e. g 0 (g%pgli)k R = \Jg^pg\ik R . 

By considering the definition of (j) 1 ^ we find 


5»PSa Te io l °*^ A ^ dt vV\ \/pTe J ° l ' / p Acdt ) = (g^i^g^ 1 ^g%pglik R ) 


,A" dt’ r -* ™ - /* i ^A^dt’ 


inducing that 


Te = \Jg<apg£g<&yfp 

I Te - i sm’Vp— R > dt ' = k R 


(277) 


(278) 


It follows from the last equality that A R ^(g 1 p) = A R (gpg t) •<=>• A^(g,p) + A R (p) = A R (gpg i). We have clearly 
A R (p) = A R (p), the “object” right connection is the Uhlmann connection. From the first equality, we see that the 


gauge potential of the “arrow” right connection is such that 

A R (g,p) = d (Vgpg^gVp) (\ / gpg ] pVp) (279) 

= dy/gpgi' \Jgpg^ + sjgpg 1 dgg~ x y/gpgi + \[gptf g^*^.{p)g~ Vs'PS' 1 ' (280) 

= Vgpg ] (^-' /: ' < A(gpg' 1 ) + dgg~ 1 g^*^(p)g~ 1 ' S j v gpg^ (281) 

T T 

/■•^(gpgf) = g'A’Hg- 1 + dgg~ x + g'Ag{g 1 p)g~ 1 by virtue of property [3] It follows that 

A R (g,p) = -\/gpg ]k g' /: v(g,p)g~ 1 Vgpg* (282) 

The connections define the curving^ of as 

B r = dA R +A R A A R + [A r , A r ] g U 2 (Morphs, f) (283) 

and the fake curvature of as 

F r = dA R + A r A A r + B r G fl 2 (Morpln// 1 , f ) (284) 


F R (\-u A , p) is just the curvature of the Uhlmann connection, it is a measure of the holonomy with respect to the 
Uhlmann’s parallelism (relative phase factor—) for infinitesimal loops starting from p. The Uhlmann’s curvature is 
related to the quantum Fisher information matrix!--, and it measures the difficulty to realize an estimation of the 
mixed state by measurements in the nieghbourhood of p. The curving is just a correction associated with the shift 
between the two purifications of gpg\ i.e. g^fp and \/gpg 1. 


We can note that it is possible to choose another “object” right connection as the another Uhlmann like connections. 
In particular, we can choose the Sjoqvist-Andersson connection by replacing A R by A BR , in this case k R is the 
Sjoqvist geometric phase, and the rest of the discussion is totally similar. The fake curvature is then written F BR = 
dA BR + A BR A A BR + B r . The possibility of changing the “object” right connection follows from the definition of the 

parallel transport of the arrow (<?», (fo 1 \fg^pgi\k R ). If we change the definition of the right geometric phase k R (by 
passing, for example, from the Uhlmann to the Sjoqvist geometric phase), the arrows (in MorphiF and Morph.#*) 
change in consequence, since ga is generated by + y/pk R k R ~ l ^fp*. 

Table fVD summarizes the different data of the connective structures. 
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TABLE VI. The different data of the connective structures of S? 1 and 


differential form 

symb. 

side 

degree 

base space 

values 

interpretation 

left gauge change of the 1st kind 

9 

L 

0 

£*(' Hs) 

G 

unitary and SLOCC operations on S 

right gauge change 

k 

R 

0 

T>\Us) 

K 

unobservable reconfigurations of A 

left gauge change of the 2nd kind 

Vk 

L 

1 

E'CHs) 

¥ 

unobservable reconfigurations of A 

left gauge potential 

a 

L 

1 

£ 4 C Hs) 

0 

generator of the operator valued geometric phase 

potential-connection 

t 

*21 

L 

1 

V\Hs) 

0 

generator of the operator valued geometric phase 

right object gauge potential 

A? 

R 

1 

V i {H s ) 

r 

generator of the Uhlmann/Sjoqvist geometric phase 

right arrow gauge potential 

A r f 

R 

1 

G x V^Hs) 

r 

unitary and SLOCC operations on S 

left curving 

B l 

L 

2 

£*( Hs ) 

¥ 

measure of the kinematic decoherence 

left fake curvature 

F l 

L 

2 

T.\U S ) 

Sc 

measure of the transitions outsite Sc(^) 

two-sided left fake curvature 

pBL 

L 

2 

£*( Hs ) 

¥ 

measure of the non-invariance 

right curving 

b r 

R 

2 

G x V\H S ) 

r 

correction needed by the arrow structure 

right fake curvature 

pR 

R 

2 

G x V\Hs) 

t 

measure of the quantum estimation difficulty 

two-sided right fake curvature 

pBR 

R 

2 

G x V\Hs) 

V 

measure of the non-invariance 


V. CONCLUSION 


The Lindblad equation describes a quantum system S in contact with a very large environment (a reservoir) 1Z. 
From a theoretical point of view, if d> £ Hs (g> H-r is the state of the system plus the reservoir, p = tvun |<I>))(( , I > | 
obeys (under some assumptions^) to the Lindblad equation. But in practice, $ is unknown because of the very large 
number of degrees of freedom of 7Z (the partial trace models the forgotting of the “informations” concerning 7Z). 
The purification process permitts to introduce 'F € Hs ® Ha suc h that p = tr- H ^|\F))((\F|, where the ancilla A plays 
the role of a small effective environment. The price to pay is the nonlinearity of the Schrodinger equation governing 
'F. The purified dynamics presented in this paper can permitt to study the dynamics of open quantum systems 
submitted to relaxation processes, with the tools used for pure states. In particular we have shown how the geometric 
phases appear in this context, unifying some concepts previously introduced by different authors. The geometric 
structure involved by the purification and the geometric phases, is richer than the case of the pure states, and needs 
the use of the category theory. We hope that these facts can help to understand dynamics of open quantum systems 
and enlighten the position of the density matrix theory in the landscape of the geometric and gauge structures in 
theoretical physics. The role of the category theory in quantum physics is an intriguing question, and it could be 
interesting to study the relations of the structure presented in this paper with the categorical structures proposed in 
other quantum problems 4 ?-—. 
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